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Abstract 

The  original  motivation  for  this  study  comes  from  the 
climatological  observation  of  three  semipermanent   sub- 
tropical highs  in  each  of  the  Northern  and  Southern  hemispheres 
and  its  explanation  in  the  work  of  H.  J.  Stewart,  which  was 
extended  by  G.  K.  Morikawa,  E.  Swenson  and  A.  S.  Peters. 

A  geostrophic  approximation  has  been  formulated  by  A.  S. 
Peters,  for  the  one  layer,  hydrostatic  atmosphere,  over  a 
rotating  sphere.   In  this  geostrophic  model,  the  coriolis 
parameter  is  treated  as  a  constant;  the  height  of  the  layer 
and  its  tangential  velocity  components  are  expressed  in  terms 
of  a  stream  function. 

Singular  solutions  composed  of  a  finite  number  of  vortices 
that  are  in  equilibrium,  may  be  found.  The  stability  under 
perturbation  of  the  following  four  types  of  configurations 
of  "geostrophic"  vortices  isstudied: 

(a)  N  vortices  of  equal  strength  m  equally  spaced  on  a  circle 
of  colatitude  9; 

(b)  N  vortices  of  equal  strength  y  equally  space  on  a  circle 
of  colatitude  6,  one  vortex  witli  strength  y.  at  the  pole; 

(c)  the  same  distribution  as  in  (a)  but  satisfying  the  boundary 
condition:  vanishing  normal  velocity  component  at  the 
equator; 

(d)  the  same  distribution  as  (b)  but  satisfying  the  boundary 
condition:  vanishing  normal  velocity  component  at  the 
equator. 
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In  this  study  we  can  distinguish  two  main  parts: 

(i)   The  solution  of  an  algebraic  eigenvalue  problem  that 
arises  from  a  linearized  stability  analysis  of  the 
nonlinear  kinematic  equations  of  vortex  motion; 

(ii)  The  numerical  solution  of  the  nonlinear  equations 

of  vortex  motion,  which  confirms,  in  some  cases,  the 
validity  of  the  linearized  analysis  and  yields  the 
magnitude  of  the  perturbed  motions. 
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0 .   Introduction 

The  study  of  geostrophic  vortices  was  begun  by 
H.  J.  Stewart  [1]  in  an  attempt  to  explain  the  occurrence 
of  the  three  so-called  semi -permanent   high  pressure  cells, 
in  each  of  the  Northern  and  Southern  hemispheres.   In  this 
work,  Stewart  represented  the  earth  by  a  rotating  tangent 
plane  and  the  large  scale  closed  isobaric  systems  of  the 
atmospheric  layer  by  discrete  rectilinear  vortices,  which 
are  defined  by  the  Bessel  function  Kq .   Later  on,  G.  K. 
Morikawa  and  E.  V.  Swenson  [2]   studied  in  greater  detail 
the  stability  of  these  single  layer  rectilinear  geostrophic 
vortices  defined  on  a  rotating  plane  with  constant  Coriolis 
parameter.   The  rectilinear  vortices  simulated  the  large 
scale  motions  of  the  atmosphere  over  the  earth  by  means  of 
projection  from  a  rotating  tangent  plane. 

A.  S.  Peters  [3]  found,  upon  neglecting  the  variation 
of  the  Coriolis  force  with  colatitude,  that  a  geostrophic 
vortex  on  a  rotating  sphere  is  given  by  a  singular  spherical 
harmonic  of  degree  v  and  order  0.   A  brief  description  of 
this  derivation  [3]  follows: 

Let   S   be  a  sphere  of  radius   a,   that  rotates  with 
constant  angular  velocity  oa   about  a  polar  axis.   Let  p  be 
the  distance  from  any  point  in  space  to  the  center  of  the 
sphere   S   and  6 ,  (})   be  the  colatitude  and  longitude  of  a 
point  on  the  surface  of  S. 
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Let   p  =  a   and   p  =  a  +  h(()),9,t)   represent  the  lower 
and  upper  surfaces  of  an  incompressible,  inviscid  fluid 
layer  that  is  gravitationally  attracted  by  S  and  suppose 
that  h((t),0,t)  <<  a. 

At  time   t  =  0,   the  constant  rotary  motion  of  the 
fluid   is  disturbed  by  the  creation  of  concentrated  vortices 
with  axes  that  are  normal  to  the  surface  of  S. 

The  velocity  components  of  a  fluid  particle,  relative 
to  S,  are  defined  by 

u  =  p  sin  9  -rr   (longitudinal  component  to  the  East) 

J  ft 
V  =  -p  -TT-       (latitudinal  component  to  the  North) 

w  =  -T^  (radial  component)  . 

The  basic  Eulerian  hydrodynamical  equations  defining 
the  motion  consist  of  a  continuity  equation  and  three 
momentum  equations. 

The  shallow  water,  or  long  wave,  model  is  derived  from 
the  Eulerian  equations  by  assuming  that: 

(a)  The  only  body  force  acting  is  that  due  to  the 
gravitational  potential   G   of   S   and  furthermore 
G  =  gp;   g  =  gravitational  constant; 

(b)  The  radial  momentum  equation  can  be  replaced  by  the 
hydrostatic  law,  for  the  pressure   p,   in  terms  of 
the  density   a   per  unit  surface  area,  that  is 

p((|),e,p,t)  =a«g'(h+a-p)  , 
which  implies 
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p  =  0   at  the  free  surface,   p  =  h+a; 

(c)  the  motion  is  such  that  the  nonlinear  terms  in  the 
momentum  equations  can  be  neglected; 

(d)  the  radial  velocity  is  negligible  and  the  radial 
variation  of  u  and  v  can  be  ignored. 

Under  these  considerations,  an  approximation  to  the 
motion  of  the  atmospheric  layer  on  the  surface  of  S  is 
determined  by  a  continuity  equation  and  two  momentum 
equations   [ 3] : 


(0.1)         n^  = 


t   a  sin 


9 (v  sin  6)    3u 

86      ~  9(f) 


/n  o\         9u    T       Q        g     9h 

(0.2)  TTX  -  2C0V  cos  6  =  -  ? ;r-  -—-     , 

9t  asine9(^ 


(0.3)  ^17  +  2a)U  cos  0  =  ^  tS"  / 

0 1.  a  0  o 

where 


h(4),e,t)-hQ 


,    h»  =h(c)),9,0)  =  constant  . 


ho      '    "0 

A.  S.  Peters  [3]  has  shown  that  by  neglecting  the 
variation  of  the  coriolis  parameter,   f  =  20)  cos  0   with 
colatitude,  equations  (0.1),  (0.2),  (0.3)  lead  to 

2       2  2 

/n    yi\  1    9    -^911,    1    9  Ti     a   r /I  2   ,  9  n  1 

(0-^)  ^i^   9?  "^"  Q  99  ^  —2-  -J   -   ihT  ^^V    ^   7^2    ^ 

sin99(J)     ^0  9t 

2   ^1 
=  2a   -r-^  C((}),9,0)  , 

where  w,   =  w  'cos  G,  =  f/2  =  oo  •  c   =  constant  and 
1  1  p 
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;-/^  A  «\       1     r  9(u  sin  9)  .  9v  ,   .   ^,     j-  , 
^(*'^'°)  =  a  sin  6  f  86 ""  3*  ^   ^^  ^^^  ^^^^^^ 

component  of  vorticity  at  t  =  0. 

For  a  vortex  of  constant  strength   y   concentrated  at 
((j).,0.)  when  t  =  0,  we  take 

2a^a)^  ^((J),e,0)    2a^a)^    6((l)-(j).)  6(6-0.) 

and   6   is  the  Dirac  delta  function. 

A  steady  state  solution  of  (0.4)  corresponding  to  a 

concentrated  vortex  is  such  that 

2oj,  ^{<p,Q) 
n(((),9)  =  -  -" 


gho 


where  ({;(({),  6)   must  satisfy 


2  2 

2     4a  u^  ijj 


(0.5)  —. A  Tq  sm  0  VH"  +  — o 9  ~  — :7i; — 

sm  0  86        86   gi^^Q  8(j)2      ghg 

5((|)-(}).)  6(6-6.) 


sin  6 


tp   does  not  depend  on  time  and  in  fact  is  the  stream  function; 
so  that  the  associated  steady  velocity  components  are 


^  =  a  39 


v  = 


1  dj^ 


a  sin  9  3(}) 

Equation  (0.5)  describes  the  geostrophic  motion  of  tha 
atmospheric  layer  over  a  sphere. 

It  can  be  verified  by  direct  substitution  that  (0.5) 
is  satisfied  by 
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(0.6)  'l^(<}>,e,<^.  ,e.)    =   ^_.^^L_  p^(-cos    a^)     , 


where   v    is    complex. 


^  =  -2-  2  y-ih^  - 1 '  -  =  ^-1  ' 

a   =  —  ,      d,  =  geodesic  distance  from  the  center 
of  the  vortex  to  any  point  ((}),e)  on 
the  surface  of  S. 
P  (x)  =  spherical  harmonic  of  order  v   [3] . 

A  superposition  of  vortices,  with  centers  at 
('I';,  (t)  ,  6,  (t)  )  satisfies  the  system  of  ordinary  differential 
equations  (1.3)  and  (1.4). 

The  stability  of  four  configurations  of  these  geostrophic 
vortices  on  the  sphere  is  studied  in  this  work: 

(a)  N  vortices  of  equal  strength  y  equally  spaced  on  a 
circle  of  colatitude  6; 

(b)  N  vortices  of  equal  strengh  y  equally  spaced  on  a  circle 
of  colatitude  6,  one  vortex  with  strength  Vq   at  the  pole; 

(c)  the  same  distribution  as  in  (a)  but  satisfying  the 
boundary  condition:  vanishing  normal  velocity  component 
at  the  equator; 

(d)  the  same  distribution  as  (b)  but  satisfying  the  boundary 
condition:  vanishing  normal  velocity  at  the  equator. 
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For  the  problems  under  consideration  there  are  four 
parameters: 

1.  Initial  colatitude  of  the  belt  of  vortices  (6), 

2.  Number  of  vortices  on  belt  (N) 

f 


3.  Coriolis  parameter   c 

4.  Strengths  y  and  y,. . 


■p    ~    2w 


The    linearized   stability  of   an   equilibriiom   configuration 
is    studied   for  N>2,      5°l6f^90°,       .2£c      <1.0    and 
-°°   <   IJq/v    <   °°'      Zones    of    linear   stability    in    the    (c    ,0)    plane 


and   the    (y./y,    9)    plane    are    given   in   the  Appendix. 


Numerical   integration   of   the   nonlinear  equations   of 
vortex  motion    (1.3)    and    (1.4)    is    used   to   test   the   validity 
of   the    linear   analysis,    i.e.,    to   see   when   nonlinear   effects 
are    important,    and  to   determine    the    linear   and  nonlinear 
motions.       Results   on   the   nonlinear   initial   value   problem 
are   given   in   Section    3   and   in    the   Appendix. 

These    computations    show   that   nonlinear   effects    are 
important  near   the   boundary   of    the    linear   stability    zone 
and  that   the   nonlinear   effects    are   more   pronounced  when    a 
polar  vortex   is    considered. 
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1.  Equations  of  Motion  for  Four  Configurations  of  Vortices 

In  general,  the  position  of  each  vortex  on  the  surface 
of   S   is  described  by  its  colatitude   6.   and  longitude  4).. 

The  stream  function   iK({),6)   is  related  to  the  velocity 
components   u,v   (longitudinal,  latitudinal)  by 

(1-1)   ("'-^  =   ilU  '   r^ITT  11)  =  (-  --  '^'   -^) 

The  stream  function  generated  by  the  superposition  of  N+1 
vortices  must  satisfy 

2    2 

1  ■i  '^r  1  '^2,  4a     CO,  til 

ir-pr  sin    9    iri:  +    ~—  — ir  -   r-^^— 


sin    e    86    -^"^    "    89    ■        .    2  2  gh. 

sm    Q    d<p  ^    0 

N  6((J)-()).)     6(8-6.) 


i=o     -  =i"  "i 

As  we  can  see  from  (0.6),  the  above  equation  is  satisfied  by 

(1.2)     i|j(4),e)  =  -j —. y   y.  P^,(-cos  a.) 

4a  sm  VTT   ^    i   v'      i' 

i=0 

with   cos  a.  =  cos  8  cos  9.  +  sin  8  sin  6.  cos  {6-6.)     , 
where   (4),  6)   is  an  arbitrary  point  on  the  sufrace  of  S 
and  (4).,  8.)  are  the  coordinates  of  the  center  of  the  ith 
vortex,   0  £  i  <_  N;   the  constant   y.  is  called  the  strength 
of  the  ith  vortex. 

It  can  be  shown  [3]  that  the  nonlinear  kinematic  equations 
of  motion  of  the  center  of  the  kth  vortex  can  be  written,  by 
using  (1.1)  and  (1.2),  as  if  the  vortex  particle  ('I'v'Qi,)  moves 
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under  the  influence  of  the  remaining  vortices: 

1     ^ 

(1.3)  a  sin  e^^d)^  =  4a  sin  vtt  .£q  ^i  ^v^'^^^  ^k^ 

•(sin  9,  cos  9. -sin  9.  cos  9,  cos  ((}),-(}).))  / 

1        '^ 

(1.4)  -a9,  =  -; ; y   y .  P  (-cos  a .,  )  sin  9  .  -sin  (ct,  -*  •  )  , 

k    4a  sm  vtt   .^-,  "^i   v'       ik       i      ^k  ^i 
1=0 

with 

(1.5)  cos  a.,  =  cos  9,  cos  9.  +  sin  9,  sin  6.   cos  {<i>T,-<t>-) 

IK  K         1  K         1  K    1 

for  k  =  0,1,... ,N,  and 


1.1.   Cases  (a,b) .   N  Vortices  on  a  Circle  of  Colatitude, 
with  and  without  a  Polar  Vortex. 

Consider  the  configuration  that  consists  initially  of 

one  polar  vortex  of  strength   Mq  and  N  vortices  of  equal 

strength  y,  equally  spaced  on  a  circle  of  colatitude  y, 

that  is  let 

9j^  =  Y  f    \  =  (k-1)  -^   ,         \i^=   \i,         (k  =  1,...,N) 
(1.6) 

9q  =  0  ,    't'o  =  arbitrary,    Mq   =   arbitrary. 

Inserting  (1.6)  into  (1.3),  (1.4),  (1.5)  we  find 
that  the  latitudinal  component  of  the  velocity  Vq   at 
the  pole  is 
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1  r     < 

^"  "  ~ sin  VTT    ^    .^   P^(-cos  y)  sin  y    sin  (())q-(})^)  =  0  , 


1=1 

since 

N 

I      sin  ((l)^-*^)  =  0  • 
k=l       "  ^ 

N 
The  longitudinal  component   u^   also  vanishes  since  I   cos  <^.-    0 

i=0 
For  the  circle  vortices,  we  find  that 

1        '^ 

V,  =  -z : y   y   P  (p,  .  )  sin  y  sin  w,  .  =  0 

k    4a  sin  vrr  "^  ,  ^ .,   v  ^ki       '       ki 

k=l 

by  symmetry,  and 

1  f      •  ,  V     . 

u,  =  a  sm  Y  ti   =  ";; = "it^nP  (-cos  y)  sm  y 

Tc  '  ^k     4a  sm  vtt  (_  0  v 

N 


i=l 
i^k 


where 


p,  .  =  -cos  Y  -  sm  Y  cos  w,  . 
'^ki         '        '       ki 

^ki  =  ^^-^^    -N 
Thus,  the  vortices  remain  on  the  circle  of  colatitude 
Y  but  they  move  around  the  circle   with   constant  angular 
velocity 

(1.7)   ^^   =    Q    =    :^ |mq  P^(-coS  Y) 

4a  sin  vir  '- 

N 

+ 


y   Z   ^v^Pki^  ^°^  ^  (l-cos  Wj^^)  I 


iT^k 
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1.2.   Cases  (c,d) .   N  Vortices  on  a  Circle  of  Colatitude, 

with  and  without  a  Polar  Vortex;  Satisfying  a  Boundary 
Condition  at  the  Equator. 

By  means  of  reflection  across  the  equator  we  can  look  for 
solutions  that  satisfy  some  prescribed  boundary  condition  at 
e  =  90°. 

The  boundary  condition  treated  here  is 


91 
8(f) 


=  0 

e=Tr/2 


i.e.      vanishing   latitudinal  velocity    component   at 
the   equator,       6   =   Tr/2 .  The   stream   function   for   the    case 

of   a  polar  vortex   along   with   N  other  vortices    in   the   northern 
hemisphere    satisfying   the   boundary   condition   is    given  by 

^1-^)       ^(*'^)    =    4a   sin   VTT      J,    ^i  [^v  (-^°-    Y^) -P,  (cos   y  ! )  ) 

1=0 

cos  Y.   =   COS  6  cos  0.  +  sin  6  sin  0.  cos  (d)-(|).) 

I 

cos  Y-   =   cos  0  cos  0.  -  sin  0  sin  0.  cos  ((|)-(j).)  . 

This  sum  can  be  interpreted  as  representing  a  superposition 
of  N+1  symmetrically  placed  vortices  in  the  southern  hemi- 
sphere with  strengths  that  are  the  negatives  of  those  of 
the  corresponding  vortices  in  the  northern  hemisphere. 

Equations  (1.1)  and  (1.8)  yield  the  nonlinear  equations 
of  motion  of  the  center  of  the  kth  vortex. 
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(1.9)  a  sin  9  4),  = 


N 


=  -;; ■■ •{  y  y  .  fp'  (-COS  a.,  )+p^  (cos  a!,  )}sin  6,  cos  9. 

4a  sin  VTT  i-r.Q  ij_^  ^       ik    ^      ik  -"      k      i 

-    (P^(-cos  a^j^)-P^(cos  a^j^))cos  Q^   sin  9^  cos  (({)j^-(|)^) 


+  2y,  P  (cos  29,  )  sin  9,  cos  6, 


(^•1°)  --\=    4i-T 


^ I   (P'(- 

in  VTT   .  ''-     V 

1=0  '- 


cos  a  .,  )  -P,,  (cos  a.,  ) 
ik    v'      ik 


iT^k 


sin  9  .  sin  (4),  -<i> .) 
1       k   1 


(1.11)   cos  a.,  =  cos  9,  cos  6.  -  sin  9,  sin  6.  cos  {<^,-(^.) 

IK  K         1  K        1  K   1 


For  the  case  that  the  N  nonpolar  vortices  have  equal  strengths 
M  and  initially  are  on  the  line  of  colatitude  y,     (1.5)  again 
yields   v,  =0,   l;^k<^N,   but  now 


(1.12)   4,j^  =  n  =  ^ 

4a  sin 


|MQ(Py(-cos  Y)+Pv(cos  y))  I 

in   VTT^^  -' 

+  y  cos  Y  I  I   P^(Pj^^)  (1-cos  w,^^)+  J  P^  (p^^^)  (1+cos  w^^)  j 

X —  X  iC —  X 


iT^k 


2      .  2 

Pj^j_  =  -cos  Y  -  sin  Y  cos  w^^  , 


Pi  .  =   cos  Y  -  sin  Y  cos  w,  .  , 
^ki         '        '       ki 

,,    .  .   27T 

^ki  =     ^^-^^     -N    • 


That  is,  the  vortices  remain  on  the  circle  of  colatitude 

Y  and  they  traverse  it  with  the  constant  angular  velocity  Q 

given  by  (1 . 12) . 
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1 . 3    Invariants 

The  kinematic  equations  of  motion  can  be  written  in 
a  more  systematic  way,  by  using  the  Kirchhoff  function  W 
for  geostrophic  vortices  on  a  sphere. 

For  a  configuration  of  (N+1)  vortices  and  no  boundary 
condition  at  6  =  it/2,  the  Kirchhoff  function  is  defined  by 


(1.13)   W  =  7i — —. y    y.li.  P  (-cos  0.  CO 

^     '       8a  sin  VTT  .  ^  -  ^    J      v       i 

-sin  9.  sin  0.  cos  ((f)^"!!)-)) 

1  J  X    J 

Then,  with  u,  =  a  sin  Qj^^*]^  and  v,  =  -a0,  ,  it  follows  that 


Since  W  depends  only  on  the  relative  distances  between 
the  vortices  and  is  not  an  explicit  function  of  time,  it 
is  easy  to  verify  that 

(i)   — -  =  0   and  therefore  W  is  itself  an  invariant,  i.e. 
dt 

W  =  constant; 

(ii)  and  from 

-  u,  a  0,   =  —■ TT-  -^^T-   I      sin  0,  7!  0 

^k    k      sm  0],  ^(f'k         ^ 

we  have 

a  I      y,  e,  sm  0   =   -  I      J^  =   ^    ' 
k=0   ^  ^      ^       k=0  '^'''k 

as  an  easy  consequence  of  the  way  that  the  angles   (J)j^ 

appear  in  W  in  (1.13).   Therefore  by  integrating,   we  get 

another  invariant  defined  by 
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N 
(1.14)   f ,  (0,  ,  62  ,  .  .  .  ,  9.  )  =   I   Ui^  cos  6,  =  constant. 

k=0 

Clearly   f,   is  an  invariant  of  the  motion. 

(iii)   Multiplying  both  sides  of 

^  =  a  sin  e^<^^      by     ^^^  q^      ,      sin  9^^  ^  0  ,  and 

summing  over  k  we  define 

N   g     1-cos  9 
^^•^^^   ^2(^1 \,    ^^,....^^)    =  Jg  ye^  sin  Q^ 

N 

_  • 

=  I     y,  (1  -  cos  e,  )  (j)  . 
k=0  '^  ^ 

For  the  case  with  boundary  condition  we  can  derive 
analogous  expressions  for  the  functions  W,  f    and  f _ . 

We  shall  see  that  although  the  function 
f2  (9,  ,  .  .  .  ,  0,  ,  (|)  w  .  .  .  / 1})!.)   given  by  (1.15)  is  not  an 
invariant,  its  linearization  with  respect  to  the  equili- 
brium solution  (2,1),  (2.2),  is  an  invariant  for  the 
linearized  motion. 

The  equations  of  motion  state  that  the  velocity  of 
the  center  of  any  vortex,  say  the  kth,  is  the  sum  of  the 
field  velocities  of  the  remaining  ones  evaluated  at  the 
position  of  the  kth  vortex. 

If  the  strength  of  a  vortex  is  positive  the  velocity 
field  of  the  vortex  has  counterclockwise  or  cyclonic 
rotation . 

In  this  system,  each  vortex  center  moves  like  a  massless 
particle  in  the  velocity  field  of  the  others. 
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2 .  Linear  Stability  Analysis  of  the  Nonlinear  Equations  of 

Motion 

The  linearized  equations  are  obtained  by  making  a 
perturbation  expansion  of  ('I'l^'Qi,)   about  the  equilibriiom 
solution  satisfying  the  initial  condition  (1.6). 
That  is,  for  the  kth  vortex 


(2.1) 


^k  =  ^  +  ^^k  ' 


<p^   =    (k-1)  ^  +  f^t  +  ea^^  ; 


while  for  the  polar  vortex 


(2.2) 

({)„  =  ea^   +   f^t  ; 

where   e  is  small  and   a,  and  6    are  functions  of  the  time  t. 


2.1.   Cases  (a,b) .   N+1  Vortices,  No  Boundary  Condition 
at  Equator. 

The  formulas  (2.1),  (2.2)  are  inserted  into  (1.3),  (1.4), 
(1.5).   By  expanding  in  powers  of   e   and  retaining  first 
order  terms,  we  get  the  following  linearized  equations 
of  motion: 
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2  •  '  f  ] 

(2.3)  -4a   sin   ^''^6,     =   iJq    P    (-cos   y){^    sin  w,  ,-ri    cos  w    ,  >• 

N        g  N        g 

1=1  1=1 

i?^k  i^^k 

2  .  f    ' 

(2.4)  4a    sin   vtt     sin    Yct=-yrv'iP    (-cos   y)cos    y 

k        u  I     V 


2    ]  ( 
+   P    (-cos   y)     sin   Yn?      cos  w,  ,+ri   sin  w^ 

^  ^kk^k  ^    .1^  ^i^i   -  J^  ^ki^i    '      1  1  k  1  N    , 
i^k  i?^k 

.    2    . 
(2.5)  ^^    ^^"    ^'^     (n+    ^O    =    (P^(-cos    Y)    N    cos    y 


II  2         '^  2  1 

+   P    (-cos   y)    sin   Y      I      cos  w . , )  >•    •    C 

i=l  '' 

r  N  N 

+   P    (-cos   y)       sin   Y      y      a •    sin  w.,-   cos  y      J      3 ■ cos   w. , 
V  '  '     i  '     . ^,       1  il  '     . ^T       1  il 


"  2 

-  P    ( -cos    y)    sin      Y      I      3.    cos   w.^     , 

i=l     ^  ^■'■ 

.    2    . 
(2.6)         ^a   sin   vtt       ^'^_^^)    =   _(p^(_cos   y)    N   cos   y 

"  2^2 

+    P„(-cos   y)     sin    Y       J      sin   w.,     1  ti 

V  '      .  ^ ,  1 1    '^ 

1  =  1 

,  N  N 

+   P,(-cos   y)     (sin   y      J      a ■    cos   w. ,    +   cos    y      J      3-    sin   w.,1 

V  '  I  /       V  '       .  ^T         1  il  '       .  ^T  1  il-* 

1=1  1=1 

2  ^ 

+   P,   (-cos    y)     sin   Y         13-    cos   w.,     ; 

V  '  .  ^^         1  11 
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where 


R         ' 
(2.7)   Aj^^  =  y  ^v^^ki^  ^°^  '''  ^^"  ^ki 


"  2 

+  y  ^v  ^'^ki^  ^^^  '^    '^°^  ^    ^^^   w,  .  (1-cos  w  .  )  , 


4i  =  ^  ^(Pki^  ^^"  ^  ^°^  ^ki-"  ^^v^Pki^  ^i"^Y  ^^^\i  ' 


I  2 

A,,  =  y^,  (cos  Y  P  (-cos  y)  +  sin  y  P  (-cos  y)) 


"kk    ^^0  ^^-^  '   v'      "       '   V 

N 

r     '  2  2 

+  y   2,   ^v^'^ki^  ^"^^^  "^  "*"  ^■'■^  "^  ^°^  ^ki^ 
i^k 

N     II  9         9  9 

+  y  j;   ^v^'^ki^  ^^^  "^      *^°^   ^  (1-cos  Wj^^)   , 


a 
(2,8)   B,  .   =  y  P  (p,  .)  cos  Y  sin  y  sin  w,  . 
ki     ^   V   ki       '      '       ki 


"  3 

+  y  P  (p,  .)  sin  Y  cos  Y  (1-cos  w,  .  )  sin  w  .  , 

V    K  X  KI  KX 


I  0  0 

A^^  =  -yP^(pj^^)  (sin  y  +  cos  Y  cos  Wj^^) 


"22  2 

+  y  P  (p,  ■)  sin  Y  cos  y  (1-cos  w  .  )   , 
V   ki  KX 


2 

-4a  sxn  vtt  cos  y^  + 


and 


^   =  Bq  cos  Uq  ,     n   =  6q  sin  a^ 


-16- 


2.2.   Cases  (c,d).   N+1  Vortices,  Satisfying  Boundary  Condi- 
tion  at  the  Equator. 

By  inserting  (2.1)  and  (2.2)  into  equations  (1.9), 
(1.10),  (1.11)  and  linearizing  analogously,  it  follows  that 


2         '  ' 

(2.9)  -4a  sin  vtt  3   =  y  <P  (-cos  y) -P  (cos  y) 

•|S    sin  Wj^^-n    cos   w,^i|+    I    A,^.B.+    I    B^.  (a    -a.)     ,       l<.k<N    , 
^  ^     1=1  1=1 

if^k  i^Oc 


2 

(2.10)     4a   sin   vtt   sin   ya,     =   -p^P    (-cos   y)-P    (cos   y) 

JC  U  I     V  V 


'COS    y\E,      cos   w    , +ri    sin  w   ^ 


-    Mq-|P^(-cos    y)+P^(cos    y)  [sin    y    j  ?cos   w^^j^+n    sin    w 


^    ^k^k    ^    },    ^ki^i    -     .?      ^ki-i    '  1    1  '^    1  N     , 

1=1  1=1 

i^^k  iT^k 


/T    ii\       4a   sin   vtt    ,'r,r\  f  fr^'  f  \ 

(2.11) (n  +  fiC)    =    Up^(-cos    y) 

I  II  II 

+    P     (cos    y))n    cos    y    +     [p    (-cos    y)-P    (cos    y)) 


2         N  2  ^0       '  1 

•    sin    Y       y      cos      w.  ,    +2    —  P     (1)  K    - 
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[p    (-COS   y)~P    (cos   y))    cos   y 


2    1    ^ 
+    (P    (-cos   y)    +   P    (cos   y))    sin   Y f    I      3.    cos    w. ^ 


N 
+    fP    (-cos   y)     -  P    (cos   Y))sin   y      J      a.    sin   w. ^ 


i=l 


2                      •  , 

(2.12) — -^  =    --^  (P    (-cos    y)+   P    (cos    y)JN    cos    y 


N  ^  2y, 


2                        2                 ^0  1 

+     fP    (-cos    y)     -   P    (cos    y)]    sin   y      J      sin    vi.-.+   P'  (1)  >    n 

1=1  ^ 

ft  f  II  II  9      1 

+s  (p    (-cos   y)-P    (cos   y))cos   y+    (p    (-cos   y)+P    (cos   Y)sin   y]  > 

N  ,  ,  N 

•    I      3.    sin    wn+    [P^  (-cos   y)-P,,  (cos   y))    sin  y      I      a.    cos  w 

i=l         1  l-L  V  V  ^^j_         1  IX 


where 


I       * 


(2.13)       aJ.    =    M(P^(Pki)     -    Pv(Pki>       ^°^    ^    ^i"    ^ki 


"^   ^'^^v^^ki^    ^   ^v^Pki^)    ^^"  ^   ^°^  ^    ^^"  \i 


y[p"(Pki)-Pv^PkiO    ^i'^^Y    cos   Y    sin   w^.    cos   w^^. 


B 


6      ^ 
k 


I      * 


i   =   ^(^v^Pki^-^^Pki^    ^^^   ^    ^°^    ^ki 


+    ^^K^Pki^    -   ^(Pki^)    ^^"^Y    ^i^^^ki    ' 
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A^j^   =    Mq    cos    y    |p^(-cos    y)    +   Py(cos   y)\ 


2       1"  "  12 

+   y^sin   Y    iP    {-cos   y)    ~   P    (cos   y)  \~    ^a   sin   vtt    cos   y^^ 


N         ,  2  2 

+   y      I      P^(P},j_)     (cos   Y   +    sin   Y    cos  w^^^) 
i=l 

"  2  2  2 

+    y      y      P    (p    .)     sin   Y    cos    y     (1-cos  w    .  ) 
.  _  -.        V        Kl  Kl 

N  ,       *  2  2 

+   y      [      ^v^'^ki^     ^^°^   ^    ~   ^^"  ^   "^^^  ^ki^ 
i=l 

*  2  2  2 

-    y      [      ^v^'^ki^    ^^^   '^    ^°^   ""^    (1+cos  w^^^) 
i=l 
iT^k 

2  2  '  2  2  2  2 

+   2y(cos   Y-sin   y)    P    (cos      -sin   y)-    8y    sin   y    cos   y 


"2  2 

P^ (cos   Y    -    sin   Y)     , 


,     .  .       * 


Bj^.    =    y    cos    Y    sin   y    sin   w,^.    P^(Pj^.)     -   P^  (P^i) 


"  3 

+   y   P    (p,  .)    sin   Y    cos   y    (1-cos    w    .)    sm  w    . 
V      kl  kl  kl 


II       *  3 

+    y   P    (p,  .  )    sin   Y    cos    y    (1+cos   w      )    sin   w  , 

V      kl  kl  kl 


A^^   =    -y   P^(pj^^)(sin   Y+    cos   y    cos    w^^) 


"22  2 

+    y    P    (p,   .)     sin   Y    cos    y     (1-cos    w      ) 
V      kl  kl 
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I   *        2       2 

-  y  P  (p,  .)  (sin  Y~  cos  y    cos  w,  .) 
V   ki        '       '       ki 


"   *       2      2  2 

-  y  ^v^'^ki^  ^"""^  '^   ^°^  '•'  (1+cos  Wj^^)   , 


and 


2       2 

'^ki  ^  ~  ^°^  '*'  ~  ^^"  '^  ^°^  ^ki  ' 

*  2.2 

'^ki  ~    "^^^  '''  ~  ^•^^  '*'  *^°^  ^ki  ' 


C  =   3q  cos  Oq  , 


n   =   3q  sin  Oq  , 


As  we  can  see,  the  linearized  equations  of  motion  have 
constant  coefficients  and  the  usual  exponential  stability 
analysis  is  possible  for  all  cases.   That  is,  the 
quantities   a   ,  B,  ,   ^  ,  n    can  be  sought  as  a  sum 
of  exponential  functions  of  the  form  e    ,  5,  =  1,...,2(N+1) 
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2.3.   Linear  Stability  Analysis  of  the  Four  Cases. 

Case  a.   N  Vortices  on  a  Circle  of  Colatitude, 
No  Boundary  Condition  at  the  Equator. 

If  we  set   Pq  =  0   in  equations  (2.3),  (2.4),  (2.7) 
and  eliminate  equations  (2.5),  (2.6),   which  involve  the 
coordinates  of  the  polar  vortex,  we  have  the  linearized 
equations  for  the  present  case. 

These  equations  can  be  written  in  the  form 


r  "^  1 

_ 

f 

B    ■ 

r    a    - 

=    L 

f    °' 

^  6   J 

^  c 

A    ^ 

^    3    ^ 

^    8 

where 


a   = 


and 


"l  "> 


^  "n 


6  = 


N 


A  =  (a^.)   = 


-  B 


a 
ki 


4a  sin  vtt  sin  y 


'^i 


4a  sin  vtt  -' 


^kk  =  °  ' 


(2.14) 


B 


(^ki^ 


ki 


4a  sin  vtt  sin  y 


C   =   (c,^.)   = 


ki 


A     2  . 

4a  sin  VTT 


N 


'kk 


=  -  1 
i=l 


'ki 


where  the  matrices    A,  B,  C   are  of  order  N. 

Furthermore   A,  B,  C   are  circular  matrices,  so  that 

it  can  be  shown  that  [8]   the  eigenvalues  of  L  are  given  by 


(2.15a) 


X  .  =  a.  +  (b.c.)^/^  , 
D     3    -        J    J 
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where 


(2.15b) 


-J-11%1,  ^.""""^ ' 

3    —    i,...,l 

•^^'^u^I.V-''"'^' 

N                      ^1,     1  ^   • 

2iTT/N 

00  =   e 

and   a.,b.,c.   are  the  eigenvalues   of  A,B,C. 

From  (2.14)  we  have  the  following  properties  of  a.,,b.,,c.,. 

IK        IK        iK 

a, ,    =    0      for   all   values   of  N 
If      N      is   even 

^l,N/2    +1   "    °    ' 

^12    "    "    ^IN    '•••'^l,N/2    "    ~    ^l,N/2    +2 


^12    = 

^IN     '•••  ''^l,N/2    " 

'^l,N/2    +2 

^12    = 

^IN    '••• '^l,N/2 

^l,N/2    +2     • 

If     N      is    odd 

^12    "    "    ^IN    '•■•'  ^l,[N/2]+l   "    "^l,[N/2]+2     ' 

^12    "         ^IN     '•••'  ^l,[N/2]+l    "       '^l,[N/2]+2     ' 

^12    "         ^IN    '•••'  ^l,[N/2]+l    "      ^l,[N/2]+2     • 

Therefore,    we   get   for    a.,b.,c.  the    following   expressions. 

3       3      3 
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a .  = 


2i  I      a^j^  sin  (u)^^  •  j) 
k=2 


(2.16)  b  .  =  b. 


P 

I 
k=2 


11  ■"  2   I   bik  ^os     (cj^^.j)  +  a(N)  cos  TTj  b^^jj^/2]+; 


c  . 


2  J   C3_j^(cos  (Oj^^.j  -  1)  +  a(N)(cos  ^  J  "D  c^  ^  [j^/2] +1 


k=2 


where 


^kl  =  ^^-^^    ^ 


P  = 


N      ^r       ^, 

2"   for  N  even 


N+1 
^  ~2~ 


for   N  odd 


a(N)  = 


1    for  N  even 


0    otherwise 


Using  (2.13),  (2.14)  and  (2.16)  we  get 


(2.17)    4a  sin  vtt  a.  =  -2ip  cos  y   I   P  (Pn^  ^i^  '^vi 

3  k=2   ^  ^ 


si 


in  (tio^i'j)  -  2iij  cos  y  sin  y   I   P  ^Pkl  ^  ^"""^  oj,  -,  (1-cos  oa,-,) 

k=2 


sin  (ti)j^j_"J)  • 


2   . 


^  ^  cos  Y 

(2.18)    4a   sin  vtt  b  .  =  -  4a   sin  vtt  —. ■-  Q. 

J  sm  Y 


P   ,  r       2 

+  2y   y   P  (p,  ,)  J^2S_X  (i-cos 
T,=2   ^   kl   [sm  Y 


'^kl'  ^°^  ^'^kl'  -'^  ^ 


+  sin  Y  [cos  w,  -.-cos  (wj^i'3)]> 


+  2y  sin  y    cos  y   J   P  (p,  ,  )  ( 1-cos  oj,  ,  )  (1+cos  (o),  ^  •  j) 

,>-^   V   kl  kl  kl 


k=2 


2 

cos  y 


+  a(N)  y  P,(p^/2  +i,i)  (il^  -  sin  y)  + 
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"^  ^   ^v^^n/2  +1,  1^  ^^^   '^         ^°^  'Yf  ^^  "•"  ^°^  ■^^^ 


(2.19)   4a  sin  vttc.  =  2ysin  y   I   P  (p,  , )  (cos  (oj  .  •  j ) -1)  cos 


CO. 


k=2 


kl 


+  2y  sin  y   I   ^v  ^'^kl^  ^^^   0Jj^-,_  (cos  (a3j^^- j  )  -  1) 


-  a(N)  y  P^(Pj^/2  +i   i^  sin  y  (cos  irj  -  1) 


(2.20)  Q.    = 


2y 


4a  sin  vtt 


J2  ''v^'k,!^  ''°^  ^  ^^"'^°^  ''kl) 


+  a(N)  cos  Y  P^(P^/2  +1,  1) 


A  Study  of  the  sign  of   c,,   shows  that   c, ,  £  0   for 

all  values  of  y ,  N,  k,  c   whenever   y  >  0  ,  which  implies 

c.  >  0  for  all  j,  y,  N,  c  . 
j  —  -'   '       p 

Therefore,  we  have  from  (2.15a,b)   that  if   b  .  <^  0 

for  all  j,   L  has  only  purely  imaginary  eigenvalues,  and 

if  b.  >  0   for  some  j,  then  L  has  complex  eigenvalues 

(except  when    c.  =  0),  the  second  alternative  would  mean 

exponential  instability  for  the  linearized  equations  if 

the  real  part  of  at  least  one  eigenvalue  is  positive. 

y  can  be  factored  from  a.,  b.,  c . ;  let 

J    3         3 


a.   =   yaj  , 


b  .  =  yb .  ,    c  .  =  yc . 
3  3  3  3 


then 


^£  =  Ma^  +  >^b;c; 
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therefore  u  affects  only  the  magnitude  of  A .  and  does 
not  affect  in  any  way  the  stability  properties  of  a  given 
configuration  of  vortices. 

If   y  <  0   we  have  to  reverse  the  inequalities  of  the 
last  paragraph. 

We  can  state  now  the  linear  stability  analysis  reduces 


to  finding  those  values  of   N, 


,  c    that  will  yield  a 
P 


negative  value  of  b .   when   y  >  0  or   a  positive  value 
of  b  .  when   y  <  0   for  all  values  of  j .   Such  values  of 
N,  0,  c   will  define  the  linearly  stable  range.   From 
(2.15a),  we  see  that  each  value  of  j  gives  two  eigenvalues 
of  L. 

If  j  =  N,  we  have   a  =  c  =  0,   which  means  that  0  is 
an  eigenvalue  of  L   with  multiplicity  2.  Summing  eq.  (2.3) 
and  (2.4)  with   y„  =  0  ,  we  have 


k=l  ^ 


N    N  N  N 

J,     J,     -ik^i  ^  J,  J,     ^ik(\--i) 

k=l  1=1  k=l  1=1 

i^^k  iT^k 


N  (    N 

M  ^    ^ik 
i=i^k=i    ^^ 

kT^i 


6, 


Integrating  we    get 


Also 


N 

I 
k=l 


N 


y      3r,   ~   constant 
k=l      ^ 


L  \  =  J.  .L  ^ik^^-^i)  ^  ^.  .L  ^ik\ 


N         N 

I         1 
k=l   i=l 

N         N 

,  ^  .     . ^ ,      ilk 
k=l    1=1 


N 

1 
k=l   i=l 


N 

V 


N      r  N 


k=l 


i=l 


^k= 


N 
-.       il 


i=l 


N 
k=l    ^ 
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By  an  appropriate  choice  of  displacements   3, 
we  can  obtain 

N 

(2.21)  ^  ^k  "  ° 

k=l  ^ 

and  then 

N 

(2.22)  I      a,  =  c  =constant  . 

k=l 

These  invariants  can  also  be  obtained  by  linearizing  f^ 

and  fy    respectively  defined  in  (1.14),  (1.15). 

If  we  make   j  =  1  in  (2.17),  (2.18),  (2.19)  we  have 

P         •  2 

(2.23)  4a2sin   vtt    a,    =    -2iy    cos   y      I      P^  (Pj^^)    sin   Wj^^ 

k=2 

2         P         "  2 

-   2iy    cos   Y    sin   y      I      ^v  ^'^kl^    ^^^   ajj^-^d-cos    ui^^) 

k=2 


(2.24)       4a      sin   vtt   b^^   =    -4a   sin   vtt    ^^^  ^    Q 


P  2  T 

r.         ,  ,         ,     cos    Y       •    •^ 

+     2U         y         P  '   (D,    -,  )      — : Sm         OJ,    ^ 

^   ,  ^^      v^'^kl'    sm   Y  kl 

k=2 

2         ^         "  2 

+    2y    sin   y    cos    y      I      ^v^'^kl^     ^^^    Wj^-^d-cos    u)^-^) 

k=2 

-  P 

(2.25)       4a      sin   vtt    c.    =    2y    sin    y      I      ^v^Pkl^     ^^°^    '^kl"^^     ^°^    '^kl 

k=2 

+    2y    sin^Y      I      Py^Pkl^     ^^"    Wj^-l(cos    (^j^j^-D 
k=2 


+    a(N)     2y   P;(p^/2    +1,    1^     ^^"   ^• 
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Using     (2.23)     and    (2.24)     in     (2.25)       we    get 

4a      sin   vtt    c,    =    4a      sin   vtt    ^'^^   ^    (n-ia, ) 
1  cos   Y  1 

Using  (2.23)  in  (2.24) , 

2  2        cos  Y 

4a   sin  vtt  b,  =  -4a  sin  vtt  — : -^  (f2-ia,  )  . 

1  sin  Y       1 

Hence 

-  iO,      and 
A,  =  a-,  +  {Q,-    ia^  )  i  = 

-*■      2a^-iQ 

For  j  =  N-1 ,  one  of  the  eigenvalues  is  A   ,  =   ifi,  i.e. 

the  period  Q      is  an  eigenvalue  of  L. 

For  Y  =  ■n^/2 ,  we  have  for  a  fixed   c    and  N, 

p 


a.  =  0   for  all   j  =  1, . . . ,N 

P 

^j  ^  ^^  k=2  ^^^''kl^  ^^°^  (.j^^-cos  cuj^i-j)- 

-  y  a(N)  P'(P^/.,  .-,  ,)  (1  +  cos  TTj) 

V  n/ Z    +1,  1 

P 

c.  =  2y  I      P^(P]^]_)  (cos  (JJj^-L'j  -1)  cos  uj^^ 


-  «(N)  y  P;(p^/2  -HI,  l)  ((-!)'  -  1) 


It  follows  that 


b,  =  0 

b^,  T  =  0   which  means  ^  =   0    . 
N-1 
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From  the  above  considerations,  we  can  conclude  that 
when   Y  =  ■n'/2 ,   we  have  only  zero  eigenvalues  for 
configurations  with  2  and  3  vortices. 

As  we  go  closer  to  the  pole  ^      increases  and 


lim  Q, 


if   y  >  0 
if   y  <  0 


when  Y  ^  0,  then  p,  .  ->   -1    for  all  i,k  =  1,...,N. 
As  the  argument  tends  to  -1 , 


p'(x)  ^  ^ill^LIL  2     and    p"(x)  ^    -   ^^^  ^^    ^ 


V  TT     X+1  V  TT      ,   ,  ,  ,  2 

(x+1) 

II 

Using  these  approximations  of  P' (x)  and  P  (x)  in 
the  expressions  for  a.,  b.,  c.  (2.17-19)   the  factor 
sin  VTT  cancels  and  we  can  conclude  that  for  small  colatitudes 
the  eigenvalues  do  not  depend  on  the  coriolis  parameter. 

The  numerical  results  listed  in  the  Appendix  suggest 
the  following  conjecture. 

Given  a  configuration  of  N  vortices,  there  is  a  colatitude 
Y,  that  depends  on  N,   such  that  for   6  <^  y   the  configuration 
is  linearly  stable  and  for   0  >  y   is  exponentially  unstable. 
This  particular  value  of   y   decreases  for  increasing  N. 
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Case  b.   N  Vortices  on  a  Circle  of  Colatitude,  1  Polar  Vortex, 
No  Boundary  Condition  at  the  Equator. 

Equations  (2.3),  (2.4),  (2.5),  (2.6)  give  us  the 
linearized  equations  for  the  case  of  N  vortices  on  a  circle 
of  colatitude   6   and  a  polar  vortex  with  no  boundary 
condition  at   6  =  tt/2  . 

These  equations  can  be  written  as 


(  5 
a 
6 


E-, 

D, 

1 

1 

Ai 

B-, 

F-, 

1 

1 

1 

c. 

A-, 

1 

1 

6 

r     6 

a 

=      ^P 

a 

.    8   J 

I    3 

where 


6  = 


r  ot. 


a  = 


a 


n 


3. 


n 


The  addition  of  a  polar  vortex  destroys  the 
circulant  property  discussed  before  and  the  eigenvalue 
calculation  has  to  be  carried  out  by  an  iterative 
procedure . 

We  have  to  consider  a  new  parameter,  namely,  the  polar 
strength   y^.   The  quotient   Mn/y   plays  an  important  role 
in  the  linear  stability  analysis. 

Some  of  the  general  properties  which  arf  valid  for  all 
values  of  the  parameters  are  the  following: 
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(a)  the  angular  velocity  Q,    as  given  by  (1.7)  is  an 

eigenvalue  of  L  . 

P 

(b)  L  has  a  double  zero  eigenvalue. 

P 

(c)  lim  Q    -  — j^ y   p'(0)  7^  0  . 

Y^tt/2      4a  sin  vtt 

In  the  derivation  of  the  linearized  equations  (2.3) -(2.6) 
mixed  spherical  and  polar  coordinates  were  used  to  describe 
the  motion  of  the  polar  vortex  and  the  vortices  on 
the  circle  of  colatitude.   This  was  necessary  in  order 
to  give  coordinates  to  the  polar  vortex  so  that  the  equation 
of  motion   is  nonsingular  at   6=0. 

If  we  sum  equations  (2.3)  and  (2.4)  over  k  or 
linearize  (2 . 14  )  ,  (2 . 15)  ,  we  get 

N  N   . 

y   3i^  =  constant    and    T   ex,  =  constant. 
k=l   ^  k=l  ^ 

That  is,  we  have  two  linear  invariants  of  the  motion 
but  they  are  independent  of  the  position  of  the  polar  vortex. 
From  equation  (1.7)  we  have  that  if 


y 


N 
cos  Y  I      Py(Pij^)  (1-cos  0)^^) 


0  i=2 


C 


* 


P  (-cos  y) 


then  Q  =   0.         Since  Q.      is  an  eigenvalue  of   L  ,  for  Vq/p  =  C 

we  have  a  quadruple  0  eigenvalue  for  all  values  of  N,  c  ,  9. 

If    p  >  0  ,  lir./\i   ^  C  ,    Q.   >_  0      and  the  rotation  is  in  the 

positive   (})   direction  (West  to  East).  If   y  >  0,  Vq/V'    <    C 
Q,    <    0. 
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Some    results   on    the    algebraic   eigenvalue   problem 
for  N   =    3    is   presented   in    this    section.         For   other   values 

of   N,    corresponding   results   are   given   in    the   Appendix. 

10      2 
The      calculations      were      done    using      y    =       .15    •    10         m   /sec. ; 


values   of      c      -    .2, 
P 


,    1.0;       0    =    5°  ,6°,  ...  ,90°    and 


VIq/p    <    0.         They   can   be   summarized   as    follows. 

N    =    3.         Mq/p    <    0.         For      a(c    )     <_  \1q/]s    <    0, 

max    |a(c    ) I    -2.0      we   have    stability    in    the   entire   Northern 

c  ^ 

hemisphere.      For       |y„/ij|     >    a(c    ),       the   stability   region    is 

reduced   to    a    cap      that   decreases   with    increasing    |  ij_/ij  |  . 


In   the    table   below  we   have    the    limit 


of    the      cap    for 


different  values   of   IJq/ij 


Table.      Of    9      for   N   =    3 
c 

Configuration   is   stable    for   0    <_  6   £ 


-2.0 

-2.5 

-3.0 

-4.0 

-5.0 

-6.0 

.2 

85° 

83° 

82° 

80° 

78° 

77° 

.3 

86° 

84° 

83° 

81° 

80° 

79° 

.4 

87° 

85° 

84° 

82° 

81° 

80° 

.5 

88° 

87° 

86° 

84° 

83° 

82° 

.6 

90° 

88° 

87° 

86° 

85° 

84° 

.7 

90° 

90° 

89° 

87° 

86° 

86° 

.8 

90° 

90° 

90° 

89° 

88° 

87° 

.9 

90° 

90° 

90° 

90° 

89° 

89° 

1.0 

90° 

90° 

90° 

90° 

90° 

90° 
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For  Mq/ij  =  C   we  have  linear  stability  for  all  9,  c  . 

Vq/v    >  0.    For   0  <  Ug/y  <  b{c  ),   max  b(c  )  =»  0.75,  we 

"      c      ^ 
have  linear  stability  for  0  £  9  <_  it/2.P 

If  Vq/v   increases,  a  belt  of  instability  appears,  and 
the  width  of  this  belt  increases  with  increasing  ii^/ii. 

In  the  table  given  below,  we  have  the  limits  of  the 

instability   belt,  6,  and  9   ,  unstable  for  Q-i  £  9  £  0->  ^ 

stable  otherwise,  for  some  values  of   c   and  y^/y. 

p      0 


Table 


Configuration  is  unstable  for  0   ;!  9  £  9p 


^ 

1 

0 

1.25 

1 

.5 

2.0 

2.5 

3.0 

4.0 

5.0 

^1 

'2 

^1 

^2 

^1 

'2 

^1 

'2 

^1 

'2 

^1 

'2 

^1 

^2 

^1 

'2 

.2 

0 

0 

11 

28 

9 

37 

9 

47 

8 

52 

8 

56 

7 

60 

7 

63 

.3 

0 

0 

10 

29 

9 

39 

9 

49 

8 

54 

8 

58 

7 

62 

7 

65 

.4 

0 

0 

10 

32 

9 

41 

8 

51 

8 

57 

8 

60 

7 

64 

7 

67 

.5 

0 

0 

10 

34 

9 

44 

8 

54 

8 

60 

8 

63 

7 

67 

7 

69 

.6 

0 

0 

10 

37 

9 

46 

8 

57 

8 

62 

8 

66 

7 

70 

7 

72 

.7 

14 

17 

10 

39 

9 

49 

8 

59 

8 

65 

8 

68 

7 

72 

7 

74 

.8 

13 

22 

10 

42 

9 

51 

8 

62 

8 

67 

8 

71 

7 

75 

7 

77 

.9 

12 

25 

10 

43 

9 

53 

8 

64 

8 

69 

7 

73 

7 

77 

7 

79 

1  .0 

12 

28 

10 

44 

9 

55 

8 

65 

8 

71 

7 

75 

7 

79 

7 

81 

♦Entries  in  table  are  in  degrees. 

Using  these  tables,  we  can  construct  the  curves 
of  linear  stability  in  the  plane  (9,yQ/y)  for  a  fixed  c 
See  Figure  A-70   in  the  Appendix. 
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Case  c.   N  Vortices  on  a  Circle  of  Colatitude, 
Boundary  Condition   at  the  Equator. 

Setting  Vq    =    ^      in  equations  (2.9),  (2.10),  (2.13)  and 
eliminating  equations  (2.11),  (2.12)  we  have  the  linearized 
equations  for  the  configuration  to  be  considered  now. 

The  equations  can  be  written  in  the  form 


a  ^ 


A2   B^ 
C2   A^ 


3 


r    a 


=   L, 


and  using  (2.9),  (2.10)  with   m»  =  0   we  have  that 


A2  =  (aj^.)  = 


(2.26)     B^  =  {b^^)    = 


-B 


a 
ki 


4a  sin  vtt  sin  y 

a 


-  A 


ki 


4a   sin  vtt 


'  \k=° 


A. 


ki 


4a  sin  vtt  sin  y 


^2  =  (^ki^  = 


B 


ki 


4a  sin   vtt 


N 


'kk 


=  -  I 
i=l 


'ki 


A^,B„,C„  are  circulant  matrices  and  therefore 
equations  (2.15)  also  hold  here. 

Using  (2.16),  we  can  get  equations  similar  to  (2.17)- 
(2.'20)  .   As  before,  we  can  observe  that   c.  ^  0  when  p  ^  0 
and  conclude  that  if  b.  >  0  for  some   j,  we  have  exponential 
instability   and  if   b .  <_  0   for  all   j,  we  have  a  linearly 
stable  configuration.    In  this  case  0  is  also  an  eigenvalue 
with  multiplicity  two. 
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Invariant  expressions  similar  to  (2.21)  and  (2.22)  can 
be  obtained  for  the  configuration  being   discussed  now. 

The  behavior  of   a  .  ,  b  .  ,  c .  as  y  "^  0   is  similar  to  that 

3    D    J 

of  case  a   and  we  can  conclude  also  that  given  N  there  is 
a   Y (N)   such  that  for   y  £  Y (N)   we  have  linear  stability. 

For  small  colatitude ,  the  eigenvalues  do  not  depend 
on  the  coriolis  parameter;  this  can  be  seen  by  using  the 
argioment  given  in  case  a. 

Even  though  the  above  properties  are  similar  to 
cases  a  and  b,  the  addition  of  a  boundary  condition  at 
the  equator  causes  some  fundamental  differences. 
First  of  all,  Q.      as  given  by   (1.12)   with   y_  =  0 
is  not  an  eigenvalue  of  L„ .   We  have  that 


fi  =  2 M  ^°s  Y 

4a   sin  vtt 


\l      P'^ip^.)     (1-cos  0^3^.) 
1=1 
i?^k 


N       * 

+   y   P'(p,  .)  (l+cos  0),  .  )  +  2P  (cos  2y) 
.  ^ ,   V   ki  ki       V 

1=1 

iT^k 


as   Y  "*"  ■'t/2,   the  terms  under  the  summation  sign  tend  to  0 
because  of  the  factor   cos  y*   But, 

lim  cos  Y  P  (cos  2y  )  =  -°° 
Y^fT/2 


so  that 


and 


lim   a  =  «>   if   y  >  0 
Y->-Tr/2 

lim   Q  =  -«  if   y  <  0 
Y^Tr/2 
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On    the    other   hand 

lim   Q  =   ±   °°      for      M    <    0 

Secondly,  as  we  will  see  later,  the  zones  of  stability  in 
the  plane   (c  ,0)   are  reduced  by  the  existence  of  a 
boundary  condition. 

Results  from  the  numerical  calculations  are  summarized 
in  the  Appendix. 
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Case  d.   N  Vortices  on  a  Circle  of  Colatitude,  One  Polar 
Vortex,   Boundary  Condition    at  the  Equator. 

Equations  (2.9)  -  (2.12)  are  the  linearized  equations 
of  motion  for  the  case  to  be  discussed,  i.e.  (N+1)  vortices 
satisfying  the  boundary  condition  v  =  0  at  6  =  tt/2  . 

The  equations  can  be  expressed: 


r  ^  ] 

• 
a 

= 

• 

^   3    J 

^3 

°3 

^3 

^3         ^3 
C3        A3 

a 
3 


=  L 


a 

3 


The  matrix   L    is  not  circular    and  the  calculation  of 
P 

the  eigenvalues  has  to  be  carried  out  by  one  iterative 
procedure . 

As  in  case  b,  we  have  a  new  parameter   Vq   and  the 
linear  stability  properties  depend  on  the  value  of  Vq/v- 

Some  of  the  general  properties  are  the  following. 

I 

(a)  the  angular  velocity  Q      is  not  an  eigenvalue  of  L  . 

P 

(b)  L        has    a   double    zero   eigenvalue. 

P 

(c)  lim   fi  =  «=   if   y  >  0  and   lim   fi  =  -°°   if   y  <  0 . 
Y^Tr/2  Y"^''t/2 

The  invariants  are  also  independent  of  the  position 

of  the  polar  vortex  and  mixed  coordinates   were  used  to 

describe  the  linearized  motion  of  the  polar  vortex  and 

the  vortices  on  the  circle  of  colatitude.    From  equation 

(1.12),  we  see  that  if 
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COS  y  |J^P^(Pj^.)(l-cos  -ki)\^/v(Pki)<l-^°^  -ki^ 
0  _    JT^^k 


P  {-cos  y)  +  P ' (cos  y) 

V  V 


* 

=  c   , 


then  ^   =   0      and   the  vortices  remain  stationary. 

*  * 

For   Mq/m  ^  c    and   y  /y  <  c    the  rotation  is  West 

to  East  and  East  to  West  respectively.   Since  Q. 

is  not  an  eigenvalue  of   L   ,  we  do  not  have  a  quadruple 

* 
zero  eigenvalue  for  Mq/P  =  c   as  in  case  b. 

As  in  the  previous  case  with  polar  vortex,  some  results 

are  given  for  N  =  3.   In  the  Appendix,  results  for  other 

values  of  N  are  given. 


N  =  3.  V-q/v    <    0.    For  negative  values  of  Mq/m   we  have 
a  cap  of  stability.   Given  below  are  values  of  the  limit, 
6   ,  of  this  cap.   (See  the  table.) 
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Table 
Values  of 


Configuration  is  stable  for  0  <  9  <  6  . 
^  —   —  c 


"p  \ 

-0.2 

-.5 

-1.0 

-1.5 

-2.0 

-2.5 

=c 

-3.0 

-4.0 

-6.0 

.2 

43° 

42° 

39° 

35° 

30° 

23° 

39° 

45° 

63° 

72° 

.3 

43° 

42° 

39° 

35° 

30° 

23° 

39° 

41° 

59° 

70° 

.4 

43° 

42° 

39° 

35° 

31° 

24° 

40° 

37° 

55° 

66° 

.5 

43° 

42° 

39° 

36° 

31° 

24° 

41° 

35° 

52° 

63° 

.6 

44° 

42° 

40° 

36° 

32° 

24° 

42° 

34° 

50° 

60° 

.7 

44° 

43° 

40° 

37° 

32° 

25° 

43° 

34° 

48° 

57° 

.8 

44° 

43° 

40° 

37° 

32° 

25° 

43° 

34° 

47° 

56° 

.9 

44° 

43° 

41° 

37° 

33° 

26° 

44° 

35° 

46° 

54° 

1.0 

44° 

43° 

41° 

38° 

34° 

26° 

44° 

36° 

46° 

53° 
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(b)  Uq/m   >   0 


For   0  <  p  /y  <  0.85,  we  have  a  cap  of 


stability  whose  limit  is 


In  the  table  below  we  have 


some  values  of  9^   for  Uq/u   =    0.25,  .5,  .75,  .85, 


Table 
Values  of 


Configuration  is  stable  for  0  <  9  <  0  , 

—   —  c 


0.25 

.5 

. —^ 

.75 

.85 

.2 

46° 

47° 

48° 

49° 

.3 

46° 

47° 

48° 

48° 

.4 

46° 

47° 

48° 

48° 

.5 

46° 

47° 

48° 

48° 

.6 

46° 

47° 

48° 

48° 

.7 

46° 

47° 

47° 

48° 

.8 

46° 

47° 

47° 

47° 

.9 

46° 

47° 

47° 

47° 

1.0 

46° 

47° 

47° 

47° 

For  larger  values  of  Pq/ij  ,  .85  <  Mq/p  <  1.5,  this  cap 
of  stability  is  broken  into  a  smaller  cap  and  a  belt  of 
stability. 

In  the  table  below  we  have  the  limits  of  the  smaller 
cap  and  the  belt  of  stability  for  values  of  Un/u    =1.0 
and  1.5. 
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Table. 
Configuration  is  stable  for  0  ^  6  f. 


)    and  e,  < 
c       1  — 


c 

p 

y^/M  =  1-0 

V^/V    =    1.5 

cap  limit 

limits  of  belt 
^1         '2 

cap  limit 

limits  of  belt 

.2 

49° 

- 

- 

8° 

43° 

52° 

.3 

49° 

- 

- 

8° 

45° 

51° 

.4 

49° 

- 

- 

8° 

47° 

50° 

.5 

49° 

~ 

- 

8° 

- 

- 

.6 

48° 

- 

- 

8° 

- 

- 

.7 

13° 

18° 

48° 

8° 

- 

- 

.8 

12° 

23° 

48° 

8° 

- 

- 

.9 

11° 

27° 

47° 

8° 

- 

- 

1.0 

11° 

29° 

47° 

8° 

- 

- 

The  belt  disappears  for  increasing  values  of  Hq/p  and  for 
Vn/u  =    2.0,  2.5,  we  have  again  a  cap  of  stability. 

Configuration  is  stable  for   0  <_  6  ^  6   .    Values  of  0  : 


.2 
.3 
.4 
.5 
.6 
.7 
.  8 
.9 
1.0 


Uq/m    =    2.0 


8° 

8° 

7° 

7° 

7° 

7° 

7° 
70 

7° 


Mq/p  =  2.5 


7° 
70 

70 

7° 
70 

70 

70 

70 

70 
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Increasing  ^n/^    '      the  new  cap,  is  broken  into 
another  cap  and  a  belt  of  stability.   Given  below  are 
some  values  of  the  limit  0  of  the  cap  and  Q-i  /9p  of  the 


belt  for  different  values  of  Mf^/y• 

Configuration  is  stable  for   0  <  6  <  6 
^  —   —  c 


14 


\^q/V 

3.0 

4.0 

6.0 

c 

p 

cap 

belt 

cap     belt 

cap 

belt 

G 
c 

^1 

4 

^   ei 

'2 

'c 

^1 

4 

.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 
1.0 

7° 

7° 

7° 

7° 

7° 

7 

7° 

7° 

7° 

60° 
63° 
64° 

63° 
63° 
64° 

6° 
6° 
6° 
6° 
6° 
6° 
6° 
6° 
6° 

63° 
63° 
62° 
65° 
66° 
67° 
66° 
67° 
68° 

67° 
67° 
68° 
68° 
68° 
68° 
67° 
68° 
68° 

6° 
6° 
6° 
6° 
6° 
6° 
6° 
6° 
6° 

67° 
67° 
66° 
66° 
65° 
65° 
67° 
68° 
68° 

73° 
73° 
73° 
73° 
72° 
72° 
72° 
72° 
71° 

If  we  compare  these  tables  with  those  of  case  b,  it 
is  clear  that  the  boundary  condition  reduces  the  zones  of 
exponential  stability  in  the  plane  (6,iJq/|j). 

The  curves  of  linear  stability  for  N  =  3  and  other 
values  of  N  are  given  in  Figures  A-69a,b,c   of  the  Appendix. 
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3 .  Niimerical  Solution  of  the  Nonlinear  Equations  of  Motion 

3 . 1   Introduction 

The  linear  analysis  of  Section  2,  established  certain 
zones  of  exponential  stability  and  instability  in  the 
(9,]j./lj)  plane  and   {e,c  )  plane  for  configurations  with 

u  p 

and  without  polar  vortex  respectively.   But,  for  the  entire 
range  of  parameters,  zero  is  an  eigenvalue  of  multiplicity 
two  and  in  some  cases  of  multiplicity  four. 

It  is  therefore  desirable  to  study  the  type  of 
solution  associated  with  the  multiple  zero  root,  and  its 
effect  on  the  linear  stability  properties,  before  describing 
the  results  of  the  numerical  integration  of  the  nonlinear 
equations  of  vortex  motion. 

First  let  us  analyze  configurations  without  a  polar 
vortex.  As  in  a  previous  section,  we  linearize  about  an 
equilibrium  solution: 

0,   (t)    =    Y   =    constant    , 

(3.1)  ^  2TT 

*^(t)  =  (k-1)  -.  ^t 

k    =    1,2,. ..,N. 

That  is,  a  perturbation  of  the  solution  (3.1)  is  considered 
to  have  the  form 

ej^(t)  -  Y  +  3j^(t)  , 


(3.2)  21T 

-I-  i^r  -t-  t 

k  =  1, .  .  . ,  N, 


,j^(t)  =  (k-1)  ^  +   Qt  +   aj^(t)  , 


where  a,  and  3,  are  small  at  time  t  =  0  . 
k       k 
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If   the   solution    (3.2)    remains    close    to   the   equilibrium 
solution    (3.1)     for    all    t,    that    is,    if    ct,   (t)    and    g    (t)  , 
k   -    1,...,N,    remain    small    for   all    t,    we   say   that    the 
configuration    is    stable. 

Assume    that   all    the   eigenvalues   of    the    linearized 
system    (2,14),     (2.26)    have   negative   or    zero   real   part,    i.e., 
the    configuration   is   exponentially    stable.      By  using   the 
invariants    defined   in   Section   2, 

N 

(3.3)  y      ^-.it)    =   C-.    =   constant 
k=l      ^  ^ 

N 

(3.4)  y      a,   (t)    =    c,'c„    =    constant    , 
k=l      ^  12 

we    see    that   if      c,    =    0,    then    the    solutions    a,   (t)  ,    6.   (t)  , 

k   =    1,...,N,    are   bounded   for   all    t,    that   is,    the    time 

increasing   components   of   the    form    (constant) 't   associated 

with   the    double    zero   eigenvalue      are   eliminated  by   the 

condition    c,    =    0    . 

If      c,    7^   0,    an    integration   of    (3.4)    gives 

N 

I      a,   (t)    =    c,    c      t   +    c^    , 
k=l      ^  ^      ^  "* 

so  that  if   c,  7^  0  ,  then  the  multiple  zero  root  yields 

time  increasing  component   in  the  longitudinal  direction. 

According  to  the  previous  remarks  about  stability,  such 

a  linear  growth  would  have  to  be  termed  an  algebraic 

instability. 
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N 
In  Section  2,  we  saw  that   c^  =   J   b.,  ,  where  b.,  is 

1=1 
given  by  (2.14).   By  using  (2.7),  (2.14),  and  (1.7),  it  is 

easy  to  prove  that 

i=l   ^     ^ 
Hence,  the  change  of  the  longitudinal  components  is 
proportional  to  the  variation  of  the  angular  velocity  with 
colatitude.   This  is  a  reflection  of  the  fact  that  when  the 
vortices  are  slightly  displaced  in  colatitude  initially, 
the  angular  velocity  of  the  equilibrium  solution  changes. 
Note  that  the  change  in  angular  velocity  is  of  the  order  of 

1  dy 
Furthermore,  as  we  will  see  later  in  this  section   and  in 
the  Appendix,  whenever  c,    ^   0,    the  vortices  tend  to 
oscillate  about  a  new  equilibrium  position,  having  a  slightly 
modified  angular  velocity  of  rotation. 

It  seems  reasonable  after  the  considerations  just 
discussed   to  define  stability  by  saying  that  the  configura- 
tion is  stable  if 

6^(t),  k  =  1,...,N,  remains  small  for  all  t; 

and  if    a,  (t)  ,  k  =  1,...,N,   remains  small  for  all  t; 

K, 

after  subtracting  a  term  of  the  form   (constant) • t. 
When  the  configuration  has  a  polar  vortex,  the 
expressions  (3.3),  (3.4)  are  independent  of  the  coordinates 
of  the  polar  vortex.   In  this  case,  the  effect  of  the  zero 
roots  on  the  solution  is  not  clearly  defined. 
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Nevertheless,  if  only  the  circle  vortices  are  perturbed 
initially,  then  the  nonlinear  calculations  show  that  the 
preceding  qualitative  description  of  the  motion  of  the  non- 
polar  vortices  also  applies  to  their  behavior  in  this  case. 
On  the  other  hand,   if  only  the  polar  vortex  is  displaced 
initially,  then  the  qualitative  behavior  of  the  resulting 
motion  cannot  be  so  easily  described. 

The  numerical  solutions   of  the  nonlinear  initial  value 
problem,  will  give  us  an  overall  picture  of  the  motions  in 
the  regions  of  linear  stability  and  instability  and  most 
important  of  all,  it  will  give  some  indication  of  the 
validity  of  the  linear  analysis. 

The  numerical  integrations ,  were  carried  out  for  a 
period  of  4000  hours  of  physical  time. 

Several  types  of  initial  perturbations  were  considered 
in  the  numerical  calculations: 

(a)  Displacement  of  the  polar  vottex  in  the  6  direction. 

(b)  Displacements  of  one  circle  vortex  in  the  positive 
and  in  the  negative  6  direction 

(c)  Displacements  of  the  circle  vortex  in  the  positive 
and  in  the  negative  <^    direction. 

The  disturbance  produced  by  displacing  only  the  polar 
vortex  is  more  nearly  symmetrical  than  is  a  disturbance 
produced  by  an  initial  displacement  in  the  cf)  direction  of 
a  circle  vortex.   In  fact  the  magnitude  of  the  latter 
disturbance  is  of  the  order  of  50%  larger  than  the  magnitude 
of  the  former  disturbance,  for  equal  initial  angular  displacement, 
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In  general,  it  is  possible  to  classify  the  motions 
obtained  by  niimerical  integration,  as  follows: 

1.  Stable  linear  motions.   This  type  of  motion   is 
typical  of  con figuratiorP  with  parameter  values  in  the 
linear  stable  zone.   They  are  characterized  by  the 
fact  that  each  component   (9  ,()),)   of  the  motion  is 
bounded,  at  all  times,  by  the  initial  perturbation 
i.e.,  the  amplitudes  of  the  oscillations  are 
proportional  to  the  amplitude  of  the  initial  perturba- 
tion (a  linear  scaling) . 

2.  Nonlinear  periodic  motions,  followed  by  an  unstable 
departure  away  from  the  initial  equilibriiom  configura- 
tion.  These  types  of  motions  occur  mostly  for  parameter 
values  in  the  exponentially  unstable  zone.   The  nonlinear 
periodic  motion  persists  for  a  finite  interval  of  time 

that  depends  mainly  on  y  /m  and  c  . 

0         P 

The  amplitudes  of  the  subsequent  oscillations  are  very 
large,  as  compared  with  the  initial  perturbation,  and  they 
are  not  related  in  a  linear  way  to  the  size  of  the  initial 
disturbance . 

Configurations  consisting  of  3  circle  vortices  and  one 
polar  vortex,  with  and  without  boundary  conditions  are 
treated  in  this  section.   Configurations  without  a  polar 
vortex  are  studied  in  the  Appendix. 
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3 . 2    3  Vortices  of  Equal  Strength  y ,  Equally  Spaced  on 
a  Circle  of  Colatitude  6 ,   One  Polar  Vortex  with 
Strength  y_  ,  No  Boundary  Condition  at  the  Equator. 

Here  we  analyze  the  numerical  behavior  of  one  case 
in  detail.  The  stationary  equilibrium  configuration  is 
defined  by  the  parameter  values 

e  =  60°,   c   =  0.2,    M  /y  =  -.597963...  . 
According  to  the  results  given  by  the  linear  analysis, 
this  configuration  is  in  the  exponentially  stable  region, 
and  zero  is  an  eigenvalue  of  multiplicity  four. 

In  Figures  3.1  and  3.2,  the  components  of  the  motion 
are  shown  for  a  polar  perturbation  of  the  form 
(gBq  =  5.0°,  ea      =   0.°)  and  (e3Q=  3.0°,  ea^   =    180°) 
respectively. 

In  Figure  3.1,  we  observe  a  steady  increase  of  the 
polar  components  in  time.   This  is  characteristic  of  the 
motions  that  result  from  initial  displacements  of  the  form 
eg   >_  3°,  ea  =0°.   On  the  other  hand,  Figure  3.2  shows 
that  a  stable  motion  with  recurrent  trajectories  occurs  for 
an  initial  displacement  in  the  direction  gBq  =  3°,  ea^   =    180° 
This  asymmetry  of  the  response  for  displacements  in  the 
eOq  =  0°  and  eoiq  =  180°  directions  arises  from  the  asymmetry 
of  the  displacements  relative  to  the  positions  of  the  circle 
vortices. 

Figures  3.3  and  3.4  show  the  resulting  motions  for  an 
initial  displacement  of  the  circle  vortex  at  (})  =  0°,  in  the 
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positive  and  negative  9  directions  respectively-   In  both 
cases  we  find  that  the  longitudinal  components  increase 
linearly  in  time.   For  small  initial  displacements,  say, 
-10°  5^  eB-,  £^10°,  the  change  in  the  longitudinal  components 
predicted  by  the  linear  theory  (see  (3.1))  is 

c,    l^=    cR       .  6.795  •  10"^  [r^^^] 
1  By      1  hour 

the    computed  nonlinear   change    is   observed   to  be   of    the   order  of 

eB,    •    7.5075    •    lO""*    [    ^^^    ]     . 
1  hour 

Clearly,  the  nonlinear  effects  are  small.   In  this  case, 
the  linear  theory  gives  a  good  qualitative  description  of 
the  motion  for  |e6-,  |  1  10°. 

In  Figure  3.5,  we  have  an  example  of  the  response  of 
the  configuration  to  an  initial  displacement  in  the  negative 
(J)  direction;  that  is   ea,  =  -2°.   Observe  that  there  is  no 
steady  drift  in  the  longitudinal  direction,  but  there  is 
some  evidence  of  increasing  components  at  the  pole.   The 
majority  of  the  integrations  performed  show  that  an  initial 
perturbation  of  a  circle  vortex  in  the  longitudinal  direction 
produces  a  "greater  deviation"  from  the  equilibrium  state 
than  does  any  other  perturbation  of  equal  initial  angular 
displacement. 
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List  of  Figures 

Each  of  Figures  3.1  -  3.5   consists  of  seven  sub- 
figures  labelled   a,b,c,...,g.   The  subfigures  are  plots 
of  the  following  quantities, 
a.    Stereographic  projection,  on  the  tangent  plane  at 

the  pole,  of  the  trajectories  of  the  centers  of  the 
vortices.   The  trajectories  are  plotted  in  two 
differently  amplified  scales.   One  scale  for  the  polar 
vortex,  the  other  for  the  circle  vortices.   The 
amplification  factors  may  be  found  by  observing 
Figures  b, . . . ,g. 


b. 

%(t) 

c. 

e^(t) 

d. 

e2(t) 

e. 

e3(t) 

f. 

(k-l) 

g- 

*n(t) 

^  -  (t)j^(t) ,  k  =  1,2,3 


For  example.  Figure  3.2b   represents  the  component  Qr,  ^^-^ 
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Figure  No.   Colatitude  of  the   Coriolis        Initial 

Circle  Vortices     Parameter    Displacement 


3.1 
3.2 
3.3 


60° 
60° 
60° 


.2 
.2 
.2 


eaQ=  5.0°,  cOq^  0° 
e3Q=  3.0°,  eaQ=  180' 
£6q=  0.01°, 
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ea.Q=    0.0°,  e6^  =  2° 


3.4 


60° 


.2 


£3q=  0.01°, 


eaQ=  0°,    e3^=  -2.0° 


3.5 


60' 


.2 


eBQ=  0.01°, 


eaQ=  0°,  ea^^  =  -2.0° 


Figure    3.1    a 
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Figure    3.2    a 
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Figure    3.3    a 
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Figure    3.5    a 
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3.3    3  Vortices  of  Strength  y.  Equally  spaced  on  a  Circle 
of  Co latitude  9 ,   One  Polar  Vortex  with  Strength  y^  , 
Boundary  Condition   at  the  Equator. 

Here  we  study  two  cases  in  detail,  for  the  coriolis 

parameter   c  =0.4. 
P 

The  linear  analysis  establishes  that  if  the  ratio 
V^r,/v    is  taken  so  that  the  circle  vortices  remain  stationary 
zero  is  an  eigenvalue  of  multiplicity  two.   Furthermore, 
for  0  <_   40°,  the  configuration  is  exponentially  stable; 
whereas  it  is  exponentially  unstable  for  0  >_   41°. 

That  is,  numerical  experiments  were  conducted  in  the 
two  cases: 

Case  (a)   9  =  35°,   c  =  .4,   y^/y  =  -0.7813...,  i.e., 
in  the  exponentiall  stable  range,  and 

Case  (b)   0  =  45°,   c   =  .4,  Hq/u   =    -0.83107...,  i.e., 
in  the  exponentially  unstable  range. 

Case  a.    In  Figures  3.6  and  3.7,  we  have  typical  examples 
of  the  resulting  motions  for  an  initial  polar  displacement, 
that  is,  they  respectively  plot  the  results  for  the  initial 
displacement  {e&      =  2.0°,  ea^   =    0°)  and  (gBq^  1.0°,  eaQ=    180°) 

We  found  that  the  motions  were  linearly  stable  motions 
for  initial  displacements  I^B^I  ±   5°. 

In  Figure  3 . 8  we  plot  the  results  for  an  initial  circle 
vortex  displacement. 

When  the  perturbation  is  in  the  positive  and  in  the 
negative  6  direction  (Figure  3.8),  we  note  that  in  this  case 
no  drift  occurs. 
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This  phenomenon  can  be  explained  by  the  fact  that  for 
the  configuration  studied 


1^  =  -  0.1637  •  10  "^  [-^] 
dy  sec 


and 


c,  |^=  -  e  .  5.895  •  10"^  ["^^^l  • 
1  3y  hour 

Therefore,  after  4000  hours  of  numerical  integration, 
the  linear  increase  is  too  small  to  be  noticeable. 

Figures  3.9  and  3.10  show  typical  motions  when  the 
initial  displacements  are  in  the  <^   direction.   The  ampli- 
tude of  the  perturbation  is  larger  in  the  case  of  a  boundary 
condition   than  in  the  case  of  no  boundary  condition. 

Case  (b) .   Figures  3.11  and  3.12  describe  the  behavior  of 
the  components  of  the  motion  when  the  initial  configuration 
is  in  the  exponentially  unstable  range.   Even  though  the 
amplitude  of  the  initial  displacement  is   small,  say 
cQq    or  eg,  =  0.1°,   we  find  that  the  motions  are  nearly 
stable  motions  for  about  1200  hours.   After  that,  we  have 
a  rapid  departure  away  from  the  equilibrium  position. 

As  was  observed  in  the  calculations  of  Section  3.2, 
an  initial  longitudinal  displacement  of  a  circle  vortex 
results  in  a  stronger  perturbation  than  that  which  arises 
from  an  initial  displacement  in  the  latitudinal  direction. 

It  is  interesting  to  note  that  a  configuration 
satisfying   a  boundary  condition  is  more  sensitive  to 
initial  displacements  than  is  a  configuration  that  does 
not  satisfy  a  boundary  condition. 
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List  of  Figures 

Each  of  the  Figures  3.6  -  3.12   consists  of  seven 
siobfigures  labelled  a,b....,g.   The  subfigures  are  plots 
of  the  following  quantities: 

a.  Stereographic   projection,  on  the  tangent  plane  at 
the  pole,  of  the  trajectories  of  the  center  of  the 
vortices.   The  trajectories  are  plotted  in  two  differ- 
ently amplified  scales.   One  scale  for  the  polar 
vortex,  the  other  for  the  circle  vortices.    The 
amplification  factors  may  be  found  by  observing 
Figures  b  ,  . . . ,g. 

b.  e^it) 

c.  e^(t) 

d.  e2(t) 

e.  e^{t) 

f.  (k-1)  ^  -  (j),  (t)  ,    k  =  1,2,3 


N 


g.   *o(t) 
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Figure  Colatitude  of      Coriolis         Initial 
Number  Circle  Vortices    Parameter,  ({>    Displacement 


3-6         35°  .4         eSq  =  2.0° 

3-7         35°  .4  6q=  1.0°,  a^=    180° 

3-8         35°  .4  6q=  0.01°,  a^^  0, 

e-L  =  -  1.0° 

3.9         35°  .4  3q  =  0.01°,  a  -  0.5° 

3-10        35°  .4  6q=  0.1°,  aQ=  180°, 

aj_=  -  .5° 

3-11        45°  .4  6q  =  0.1°,  a^-  0° 

3-12        45°  .4  6q  =  0.01°,  6-^  =  0.1° 
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4.   Numerical  Method 

The  linearized  eigenvalue  calculation  of  Section  2 
was  carried  out  by  using  tlie  eigensystem  subroutine  package 
EISPACK  [9].   The  routines  from  this  package,  reduce  the 
input  matrices   (L,  L   ,  L   ,  l'  of  Section  2)  to  upper 
Hessenberg  form  by  means  of  similarity  transformations  and 
then  computes  the  eigenvalues  by  applying  the  QR  method. 

The  numerical  integration  of  the  nonlinear  equations 
of  vortex  motion  was  carried  out  by  using  the  modified  Euler 
method  [5]  defined  as 

u   =  u^  +  At  f(t,u^) 
(4.1) 

^t+At  ^  ^t  ^  At  ^^(^^^t)  ^  f(t+At,u*)) 

The  equations  of  motion  of  the  center  of  the  polar 
vortex,   have  a  singularity  at   0-  =  0,  which  comes  from 
the  system  of  coordinates  used  to  describe  the  motion. 
When  the  initial  displacement  is  a  polar  displacement, 
this  singularity  presents  no  problem  since  9.  (t)  -^   0  for  all  t. 
In  the  case  of  an  initial  circle  vortex  displacement, 
the  longitudinal  component   (J)„   of  the  polar  vortex  is 
not  defined  at   t  =  At.   In  equations  (4.1) ,  the  function  f , 
corresponding  to  the  equation  defining  ^^    ,    has  a  multipli- 
cative factor  of  the  form  l/(sin  6„  ). 

Numerical  experiments   were  carried  out  with  two 
different  methods  to  overcome  the  above  difficulty: 
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(a)   At  t  =  0  ,  displace  both  the  polar  vortex  and  a  circle 
vortex,  in  such  a  way  that  the  polar  vortex  displacement  is 
equal  to  a  small  fraction  of  that  of  the  circle  vortex.   The 
effect  of  the  initial  displacement  at  the  pole,  is  negligible 
compared  with  the  effect  of  the  initial  circle  vortex 
displacement  if 

l-^l  <  0.1 


where  e„  is  the  initial  angular  displacement  at  the  pole  and 
e^      is  the  initial  angular  displacement  of  the  circle  vortex. 
In  this  way,   ^n^^^   "'"^  nonzero,  the  singularity  is  avoided 
and   the  resulting  motion   has  the  overall  characteristics 
of  a  circle  vortex  perturbation.   In  the  numerical  calcula- 
tions we  reported  in  Section  3,  values  of  leQ/e, |  were  taken 
to  be  0.01,  0.02  and  in  one  case  0.1.   From  the  equations 
of  motion  given  in  Section  1,  it  is  easy  to  see  that  the 
disturbance  created  by  the  initial  displacement  of  the  polar 
vortex  is  of  the  order  of  e^. 

(b)   If  the  velocity  vector  q   of  the  polar  vortex  is 
considered  to  be  in  the  tangent  plane   (x,y)   at  the  pole, 
then  its  components   (q  ,q  )   in  that  plane  are  given  by 


(4.2) 


r  "3 


X 


^  q. 


-sin  (}).   -  cos  4)_ 


cos  (pQ      -   sin  (})- 


r    u. 


V 


0 


when  the  polar  vortex  is  considered  to  be  at  6  =  0 ,  (()  =  c})^ 
and  to  have  spherical  velocity  components   (Uq,Vq). 
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On  the  other  hand,  at  t  =  0,  the  polar  vortex  satisfies 
6q  =  0 ,  'pQ   ~   arbitral  y.   We  can  calculate   Uq  =  a  sin  O^^'r)  / 
v„  =  -  aB^   by  using  the  equations  of  motion  of  Section  1. 
We  find  that 


N 

(4.3)  u„  =  -  -. ^ y   P  {-cos(0,  +eB,  )  )  sin  (9,  +e6,  ) 

0      4a  sin  vtt    .  -^  ^   v       k   k         k    k 

k=l 

•cos  ('l^k+^o'j^'^o) 

=  u    cos  ({Iq  +  V  sin  (])„ 

N 

(4.4)  v„  =  -; 7^ y   p' (-cos(0, +e3,  )  )  sin  (e,+e6,) 

0    4a  sin  vtt    ,  ^^   v       k   k         k   k 


k=l 


=  -u  sin  ({)_  +  V  cos  cj) 


sin  {(t>^+ea^-(pQ) 


Hence,  the  (x,y)  coordinates  of  the  center  of  the  polar 
vortex  may  be  predicated  at  At  =   t  from 

X  =  At  q^  , 

y  =  At  q^  . 

Using    (4.2),     (4.3)    and    (4.4),    we    see    that 

q      =    -V 
^x 

q      =      u    . 
If  we   map   the    coordinates      x,y      to   the    sphere   by  means 

of    a    stereographic   projection, 

2  2       2 

*  r4a^    -    (x^+y^)     ^ 

6      =    arccos    [ — 5 22 — ^ 

4a      +      x   +y 

4)      =   arctan    (3^)    / 

*  * 

we    find   predicted  values    for    9„(At)    =    9      and    (J).  (At)    =    cj)    . 
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Now  by  using  a  first  order  corrector  formula,  the 
polar  vortex  is  located  at   At   from 

<l>o^^  =  *o°  +  At  g  (e*,(t.*,e.°,(})^°) 

0Q^^  =  Qq^   +  At  h  (e*,(t)*,e^°,(j)^°)  ,   i  =  i,...,N, 

where   g,  h   can  be  found  using  the  equations  of  motion 
of  Section  1.   We  could  have  obtained  corrected  values 
for   0-(At)  and   c}).  (At)  ,       accurate  to  second  order 
by  using  an  iteration  of  the  form 

t+At    t  ,  ,^  ^/^  ,  At   u^+u   , 
u      =u   +Atf(t  +  -2,  — 2 * 

With  this  starting  method,  we  continued  the  original 

numerical  integration  sequence  (4.1)  until  Q^      become 

-2 
smaller  than  a  prescribed  value,  say   10    degrees,  at 

which  point  a  similar  mapping  was  needed  to  continue 

the  integration  until  the  polar  vortex  passed  away  from 

the  pole . 

Some  integrations  were  carried  out  by  using  both  alter- 
natives (a)  and  (b) .   The  results  show  no  significant 
differences.   Hence  method  (a)  was  used   primarily  owing 
to  its  simplicity. 

The  spherical  harmonic   P  (x)   that  defines  the    stream 

function  for  geostrophic  vortices  on  a  rotating  sphere, 

satisfies  the  Legendre  equation 

2 
(4.5)    (1-x^)  P;'(x)  -  2  X  P^(x)  -   ^^  l^-   P^(x)  =  0  , 

where 
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,   ,.,2  2      1/2 
2  ^  g  hQ      ^ 

V  =  -  2  +  1  A 


X  =  -  cos  a . , 


=  -  cos  e,  cos  9.-  sin  6,  sin  6.  cos  (({i^-cf)-)  / 

K         X  K         1  K    1 


(9    ,(J)„)       coordinates   of    t±ie    Ith   vortex   point. 
P    (x)    can   be  written    [7]  : 


P    (x)    -   F(y  -    iA,    \  +iX,    1;    Kr-] 


1    ^    1+4A^     rl-x^     ,      (1  +  4A^)  (3^  +  4X^)        1-x    2 

~  2'-?-'  22  2 

2  2     '4 

,     (1+4A^)  (3^  +  4A^)  (5^  +  4A^)  rl-x^ 3 

o2     ,2    ^2  ^    2    ^  •  •  • 

2    •  4    •  6 

=  1  +  a   +  a   +  . .  for   -1  <  x  <  3. 

From  the  series  expansion  of  P  (x) ,  we  see  that 

(  (2n-l)^  +  4A^)  rl-x^ 

^n+1   -   ^n  *       ,„  .2     ^  2  ^ 

(2n) 

1+4A^  ril2i1 

The    squares   of   the    even   and  odd   integers    can   be 

calculated   once      and   stored.      Furthermore,    for    fixed      A, 

2         2 

we  calculated  and  stored    "~— 5 .   The  calculation 

(2n)^ 
of  one  term  of  the  series  a   requires  only  2  multiplica- 
tions and  therefore  the  recursive  calculation  of   a    is 
an  efficient  way  to  calculate   P  (x)  .     Since   P^  (>^)  has 
a  logarithmic  singularity  at   x  =  -1 ,  for  values  of  the 
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argument   -  0.995  <  x  <  -0.94,   i.e.   5°  <  a.,  <  20°,  about 

—   —  'X/   ik  'X/ 

400  terms  of  this  series  were  necessary  to  get  a  relative 

I 

accuracy  of  better  than  1/10  00  for  the  value  of  P  (x) . 

The  first  derivative  of   P  (x)   was  evaluated  by  using 
the  derivative  of  the  series  for  P  (x)   while  the  second 
derivative  was  calculated  by  using  (4.5)  and  the  already 
computed  values  of  P   and  P'. 

The  numerical  integrations  were  monitored  by  means  of 

the  invariant  (1.14), 

N 

I(t)  =  I      y^  cos 

k=0 

In  general,  the  maximum  relative  change 

I I(t)-I(  0 )  I 
E  =  max  l-y^o-j 1 

-4 
was  bounded  by  .I* 10   .    For  initial  displacements  less 

-7 
than  5°,   E  was  bounded  by  0.1  •  10 
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5 .   Comparison  of  Geostrophic  Vortices  on  the  Sphere 

with  Recitlinear  Geostrophic  Vortices  on  a  Rotating  Plane. 

It  is  of  interest  to  make  a  qualitative  comparison 
between  the  rectilinear  geostrophic  vortices  studied  by 
G.  K.  Morikawa  and  E.  V.  Swenson  in  [2]  and  the  geostrophic 
vortices  on  the  sphere  studied  in  this  work. 

A  single  rectilinear  geostrophic  vortex  on  a  rotating 
plane  is  defined  by  the  stream  function  [2] 

"^0  =  2?^0(^  l^-^o'^ 
where   Y  =  strength  of  the  vortex 


K  =  2(j//ghQ  , 

u>  =  angular  velocity  of  the  plane 
! r-r„ !  =  distance  from  an  arbitrary  point  on  the  plane 

to  the  center  r„  of  the  vortex 
K_(x)  -   Bessel  function  of  order  0. 

The  stream  function  for  a  distribution  consisting  of  one 
vortex  with  strength   Yq   at  the  origin  of  the  system  of 
coordinates  and  N  vortices  of  equal  strength  y ,    equally 
spaced  on  a  circle  of  radius   r,   centered  at  the  origin, 

is  given  by 

Yn  '^ 

^  =  ^^i-    Ir-r^l)  +  ^7  .I^Ko(K  |r-r.|)   [2]. 

Since  in  [2]  the  radius  of  the  circle  is  normalized  to  1, 
for  the  plane  problem  tliere  are  three  independent  parameters, 
(1)   N,  the  number  of  vortices;  (2)  k  ,  scale  parameter, 
(3)   Ya/y^   quotient  of  strengths. 
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In  the  spherical  case,  for  a  distribution    consisting 
of  one  polar  vortex  with  strength   )j„  ,   N  vortices  of  equal 
strength   p,   equally  spaced  in  a  circle  of  colatitude   9, 
no  boundary  condition  at  the  equator,  we  have  4  independent 
parameters : 

(1)   N,  number  of  vortices,  (2)  6,  colatitude  of  the  circle 
vortices,   (3)  y/^/M/  quotient  of  strengths,  (4)  c   ,  coriolis 
parameter. 

The  first  observation  is  that  the  spherical  problem 
depends  on  one  more  parameter  and  a  direct  parametric 
equivalence  cannot  be  easily  established. 

However,  two  important  associations  between  the  stream 
functions  can  be  made: 
(1)  Peters  [3]  noted  that 

1^"^  4  sin  vu  Pv^-^^^  ^1^  =  -  47  S[<yr^+rl-2r   r^cos  ((f,-(D^) 


4 

sin 

203 

VTT 

•c 

whe 

re     K 

p 

a  =  radius  of  the  sphere. 


(2)  As  it  was  stated  in  Section  2,  for  small  values  of  a,  , 
the  stream  function  for  one  vortex  in  a  rotating  sphere 

i)   = r-= P  (-cos  a,)  , 

^    4  sm  VTT   V       1 

cos  a,  =  cos  9  cos  9,  +  sin  9  sin  Q-^   cos  ((})-(|),) 

becomes  independent  of  c  .   Furthermore, 

P 
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1-cos    a, 
P^(-cos    a^^)    '^^   A    log    ( 2 )     • 

On    the   other  hand,    for        k       close    to   0      we   have    that 

ir-r    I 
Kq(  k  |r-rQ|)    -v.   B    log   -^ , 

so  that  both  stream  functions  have  the  same  type  of 
logarithmic  behavior  for  values  of  the  argument  close 
to  their  respective  singularities. 

Since  the  two  problems  have  the  same  type  of  symmetries 
and  geometric  properties,  the  matrices  of  the  linear 
problem  have  the  same  form.  This  fact  makes  some  basic 
properties  of  the  linear  analysis  common  to  both  cases: 

(a)  zero  is  an  eigenvalue  of  multiplicity   two  and  in  some 
cases  of  multiplicity  four. 

(b)  The  equilibrium  angular  velocity  9.      is  an  eigenvalue 
of  the  linear  matrix. 

(c)  There  are  two  linear  invariants.   When  a  polar  vortex 
is  present,  they  are  independent  of  the  coordinates 
of  the  polar  vortex. 

(d)  Two  distinct  families  of  eigenvalues  with  positive  real 
part  are  evident,  one  family  is  pertinent  to  the  circle 
vortex  and  the  second  family  to  the  center  vortex. 

(e)  In  particular,  for  N  =  2,   y^/y  =0   is  in  the  unstable 
range.   If   Va/i^  ^  ^'   ^^®  configuration  is  exponentially 
unstable  for  all  values  of  the  parameters. 

(f)  For  N  7^  3   the  configurations  are  exponentially 

unstable  for  all  values  of  the  parameters  when 

V'q/V'      =    C   (Section  2)  . 
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Some  important  differences  can  be  listed  as  follows; 


Planar  Case   [2] 

Spherical  Case 

(a)  vig/y  <  0 

(a) 

Vq/v    <    0 

N  =  3  is  exponentially 

N  =  3  is  exponentially 

stable  for   all  values 

stable  for  some  values 

of  <  . 

of  c   and  6. 
P 

4  <  N  <_  6 ,  exponentially 

N  >^  4  exponentially  stab 

stable  for  some  values 

for  some  values  of  c  , 

P 

of  K 

N  >  7,  exponentially 

unstable  for  all  values 

of  K  . 

(b)  Mq/m   =   0 

(b) 

Mo/y  =  0 

N  =  2,3,4,5,  exponenti- 

N  =  2,3,  exponentially 

ally  stable  for  all 

values  of  K  . 
N  =  6,  exponentially 

stable  for  some  values 

of   K. 
N  >_  7,  exponentially 

unstable  for  all  values 

of   K.  However,  for  N=7, 

K  =  0  is  a  neutral 

stability  point. 


for  all  values  of  c  , 

ir 


N  >  4,  exponentially  stable 


for  some  values  of  c  , 

P 
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Consider  a  configuration  of  three  vortices  of  equal 
strength   y   equally  spaced  on  a  circle  of  colatitude, 


one  polar  vortex  with  strength 


y^  ,   no  boundary  condition. 


and  a  choice  of   Mq/m   so  that  the  circle  vortices  remain 
stationary. 

Bauer  and  Morikawa  [2a]  established,  by  using  a  plane 
model,  that  the  above  system  of  vortices  can  be  used  to 
represent  the  climatological  semi -permanent  highs  whenever 
1.45  <_  K  <_  5 . 0 .   These  values  of  k     imply  that  the  ratio 
q,  quotient  between  the  long  period  and  short  period  of 
the  system,  takes  on  values  between  2  and  32.85. 

In  the  spherical  case,  values  of  q  between  2.0  and 
10.0  are  found  at  colatitudes   55°  ^  6  ;^  75°,  which  is 
approximately  the  geographical  region  where  the  semi- 
permanent highs  are  observed. 

In  particular,  if  we  take  the  magnitudes  of  ij„  and  y 
in  such  a  way  that  2000  Km  away  from  the  center  of  a 
circle  vortex,  the  speeds  are  about  10  [m/sec] ,  we  obtain: 
(a)   Planar  case 


e  =  66.5° 

K      =  3.75 

q  =  12 

]Aq/u    =    -    0.081 

Short  period:  1  year 


6  =  60° 
K  =  3.2 
q  =  7.7 

-  0.121 
short  period:  0.52  years 


Mg/li 
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(b)       Spherical   Case 

9    =   66.5° 

c      =    1.0 
P 

q  =  8.45 
y^/y  =  -.119 
short  period: 
1.6  years 


e   =  60° 

c   =  1.0 
P 

q  =  6.89 
M^/V=    -    .144 
short  period; 
0.4  years. 


In  Figure  5.1  we  have  plotted  |yQ/y|  =  \c*\    versus  q 
for  the  following  cases: 

(1)  Plane  mod'^il,  curve  A.   The  left  end  of  the  curve 
corresponds  to  k=  3.6,  while  the  right  end  corresponds 
to  K  =  2.0. 

(2)  Spherical  model  with  circle  vortices  at  6=  70°,  curve  B. 
The  left  end  of  the  curve  corresponds  to  c   =  1.0  while 
the  right  end  corresponds  to  c   =  .2. 

(3)  Spherical  model  with  circle  vortices  at  0  =  60°,  curve  C. 

Left  end:  c   =  1.0,   right  end:  c_  =  .4. 
P  P 

(4)  Spherical  model  with  circle  vortices  at  6  =  50°,  curve  D. 

Left  end:   c   =  1.0,   right  end:   c   =  .5. 
p  P 
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A.   Appendix 

1.    Numerical  Results  on  the  Eigenvalue  Calculation. 

As  it  was  stated  before,  the  strength   y   does  not  have 
any  effect  on  the  linear  stability  properties.   For  the 
majority  of  numerical  experiments,  it  was  chosen  to  be 

M   =   .15  X  10     ]  . 

seg 

For  this  value,  the  speed  of  the  velocity  field  20°  away 

from  a  center  of  a  circle  vortex  is  between  8[m/sec]   and 

20 [m/sec] . 

Because  of  meteorological  considerations,  the  following 

choice  of  h^ ,  a,  g  was  made 

hg   =   8500  [m] 


a   =   6370  [km] 

~2 


g    =   9.806  [-^^] 


seg 

Computations  were  made  for  N  =  2,3,. ..,8,  c  =  .2,. 3,..., 
1.0   and   0  =  5°,6°, . .  .  ,90°. 

i)    N  Vortices  on  a  Circle  of  Colatitude,  No  Polar  Vortex, 
No  Boundary  Condition   at  the  Equator.  (2  <_  N  <_   8) 

For  N  =  2,3  we  have  exponential  stability  for  0  <,  9  ;!  y* 

For  N  >  3,   a  polar  cap  of  stability  is  characteristic  for 

each  N.   Table  A.l  shows  the  values  of  the  limit   6    of  the 

c 

stability  cap.   As  we  can  see  there,  for  fixed  N,   6   decreases 
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with   increasing   c   and  for  fixed  c  ,  9   decreases  with 

P  P    c 

increasing  N.   Table  A. 2   shows  the  values  of   6    with 

c 

c   =  cos  0,  i.e.,  using  the  true  value  of  the  coriolis 
parameter. 


Table  A. 2   Values  of  6   for  c   =  cos  6. 

c      p 

Configuration  is  stable  for  0  £  G  £  G  . 


N 

'o 

4 

50.4093° 

5 

37.9224° 

6 

17.5225° 

7 

9.9952° 

8 

8.4092° 

Let   w(N,0,c  )  =  max  | Im  A. I . 
P      i   '     i' 

This  function  will  give  us  an  overall  view  of  the  behavior 
of  the  maximum  frequency  with  respect  to  N,  G,  c  .   Its 
properties  are  the  following: 
1.    For  fixed   N   it  is  continuous  with  respect  to  G  and 


c   for  0  <  0  <  tt/2  ,   .2  <  c   <  1.0. 
P  -   P  - 

2.  lim  aj(N,G,c  )  =  °°   for  all  N,  c   . 
6^0         P  P 

3.  fi  5  co(N,0,c  )  for  all  N,G,  c  . 

P  P 

4.  1^  (N,G,c  )  <  0     for  all  N,  c    and  for  8  <  9  . 

do  p   —  p  C 


<_   9   the  eigenvalues  of  L  have  the  following 


For    _ 
distribution . 

1.  1  pair  of  0  eigenvalues; 

2.  +    i    Q,    ; 
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3.    n-2  pairs  of  imaginary  complex  conjugate  eigenvalues. 

For   0>0   +0.00  01°  we  have 
—  c 

1.  Double  zero  root 

2.  +  i  f^ 

3.  Eigenvalues  of  the  form  ±a,,±a  tib^/  a,,a2/b   real. 

In  Figure  A-1,  the  curves  of  linear  stability  with 

respect  to  9  and  c   have  been  drawn.   For  N  fixed,  we  have 

P 

linear  stability  to  the  left  of  each  curve. 

As  we  can  see,  when  N  increases,  the  curve  becomes  more 

parallel  to  the  c   axes  and  therefore  the  value  6^  does  not 
'^  p  c 

depend  so  much  on  c  .   This  result  stems  from  the  fact 
stated  earlier:       when   0  is  small,  the  eigenvalues  are 
independent  of  c  .   A  graphic  representation  of  max  | Re  A . | 
versus  colatitude  (0)  and  Im  |A.|,  i  =  1,...,N  for  different 
values  of  the  parameters  are  given. 
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ii.   N  Vortices  on  a  Circle  of  Colatitude,  No  Polar  Vortex, 
Boundary  Condition   at  9  =  tt/2  .    (2  <_  N  £  8) 

For   N  >  2,   we  have  a  cap  of  stability.   In  particular 
for  N  =  2,3   we  do  not  have  exponential  stability  for 
0  <_  0  <_  Tr/2   as  the  previous  case  but  a  polar  cap.  For  N  >  4, 
we  have  a  polar  cap  of  stability,  but  this  cap  is  smaller  than 
that  of  (i)  for  each  N. 

In  Table  A. 3  we  have  the  limit   0   of  the  cap  of 

c  '^ 

stability  for  different  values  of  N  and  c  .   In  Table  A. 4  we 

P 

have  the  values  of  0   for  c  =  cos  0. 

c      p 

If  we  compare  the  entry  for  N  =  8  in  Tables  A. 2  and  A. 4 
we  can  see  that  they  are  equal  up  to  four  decimal  digits ; 
the  same  happens  in  Tables  A.l  and  A.  3  for  N  =  8,  c   =  .9  and 
c  =1.0.    This  suggests  that  when  y   is  small,  the  stability 
properties  of  both  cases  are  similar,  independent  of  the 
existence  of  a  boundary  condition. 

We  observe  that  a.,  b.,  c,  (2.16),   depend  on 

:   J   D 

p'(-cos  Y  -  sin  Y  cos  co,  .  )  , 

V  ki 

2        2 
P"(-cos  Y  -  sin  Y  cos  co,  .  )  ; 

_  V  Kl 

and  2         7 

P'(cos   Y  -  sin  Y  cos  to  .  )  , 

V  kl 

2         2 
P"  (cos   Y  -  sin  Y  cos  co,  .  ) 

V  kl 

where  the  first  pair  of  functions  comes  from  the  vortices 

on  the  northern  hemisphere  and  the  second  pair  from  the  images 

in  the  souterh  hemisphere. 
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As   Y  ^  0,  the  first  pair  goes  to  +  °°  respectively 
and  the  second  pair  tends  to  a  constant;  hence  the  contribution 
of  the  images  is  very  small  compared  with  the  contribution 
of  the  vortices  on  the  northern  hemisphere,  and  therefore 
the  linear  stability  properties  are  the  same  for  both  cases 
when   Y   is  small.   From  this  we  can  say  that  for  large  N, 
the  cap  of  stability  tends  to  be  the  same  in  both  cases. 


Table  A.  4.    Values  of  9   for  c  ==  cos 

c      p 


Configuration  is  stable  for   0  ^ 


N 

'o 

2 

38.4354° 

3 

45.3942° 

4 

41.7936° 

5 

36.2536° 

6 

17.5151° 

7 

9.9951° 

8 

9.4092° 

9 

8.1247° 

It  is  also  worth  mentioning  that  ^   is  larger  when  we 
have  a  boundary  condition;  this  is  due  to  the  increase,  in 
the  longitudinal  component  of  the  velocity,  that  comes 
from  the  images  - 
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In  Figure  A. 26,  we  have  the  curves  of  linear  stability 

with  respect  to  9  and  c  .   For  fixed  N,  we  have  linear 

P 

stability  to  the  left  of  each  curve. 

If  we  compare  Figures  A. 1  and  A. 26,  we  see  that  for 

N  =  6,7,8,  the  curves  have  the  same  form;  if  N  <  6  there 

are  noteworthy  differences.   For  N  =  2  and  5,  they  intersect; 

also  note  that  for  N  =  2,3,   dc  /dB  >  0  and  for  d    >  3, 

dc  /d9  <  0.     For  N  =4,5,   the  curves  are  shifted  to  the 
p 

right. 

In  general  we  can  conclude  that  the  addition  of  a 
boundary  condition  makes  the  configurations  more  unstable 
in  the  sense  that  the  zones  of  stability  are  reduced. 

A  graphic  representation  of  max  | Re  A. |  versus 
colatitude  and  Im  | X . |  ,  i  =  1,...,N,  for  different 
values  of  the  parameters   are  given. 
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iii)    N  Vortices  on  a  Circle  of  Colatitude, 

One  Polar  Vortex,  No  Boundary  Condition  at  9  =  j.  (2<I'I£8) 

Calculations  were  made  for  N  =  2,3,. ..,8,   c   =  .2, ...1.0, 

P 

9  =  5°,6°,  .  .  .  ,90°   and  values  of  Mq/u   =   c*,  Mq/m  >  0, 


Mq/p  <  0. 
as  follows, 
(a)    N  -  2 


The  general  results  can  be  summarized 


If  M^  y  >  0   we  always  have  exponential  instability 

for  all  values  of  G ,  c  ,  yn/M- 

p    0 

If  Mq/m  <  0   we  have  linear  stability  in  a  belt  whose 

limits  are  6   and  tt/2  . 
c       ' 

The  table  below  qives  some  values  of  6   for 

^  c 

c   =  .2,...,  1.0   and  for  y  /y  =  -0.125,-0.25,-0.5,-1.0,-1.25, 
-1.5. 


-0.125 

-0  25 

-0.5 

-1.0 

-1.25 

-1.5 

.2 

87° 

83° 

74° 

46° 

11° 

<  5°  ' 

.3 

87° 

83° 

74° 

47° 

11° 

<  5° 

.4 

87° 

83° 

74° 

48° 

12° 

<  5° 

.5 

87° 

83° 

74° 

50° 

13° 

<  5° 

.6 

87° 

83° 

74° 

51° 

18° 

<  5° 

.7 

87° 

83° 

74° 

52° 

24° 

<  5° 

.8 

87° 

83° 

75° 

52° 

29° 

<  5° 

.9 

87° 

83° 

75° 

53° 

32° 

<  5° 

1.0 

1 

87° 

83° 

75° 

54° 

35° 

<  5° 

Table  of  9   for  N  =  2.  Configuration  is  stable  for  6  <Q<_   tt/2. 
For  \i(^/u   =   c*  the  configuration  is  unstable. 
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As  we  can  see  when  |)j-/y|  increases, 


decreases  and 


the  configuration  becomes  stable  for   q*  <  q  <  90°   0*  <  51 

(b)    N  =  3. 

Pq/m    <    0.       For    a(c    )     <_  y^/y    <    0,    max    Ia(c    )  |    ~    2.0, 
we   have    stability   in    the   whole   northern  hemisphere. 


For    |y_/y|     >    a(c    ),    the    stability    region    is    reduced   to 

a   cap    that   decreases  with   increasing    |  y^^/y  |  .       In   the    table 

below   we   have    the    limit    9      of    the    cap    for   different   values 

c 

of  y^i/y  and  c  . 
0        p 

Table  of  0   for  N  =  3.  Configuration  is  stable  for  0  <  6  <  6  . 
c  ^  —   —  c 


■2.0 


-2.5 


-3.0 


-4.0 


-5.0 


-6.0 


.2 

85° 

.3 

86° 

.4 

87° 

.5 

88° 

.6 

90° 

.7 

90° 

.8 

90° 

.9 

90° 

1.0 

90° 

83° 
84° 
85° 
87° 


82° 
83° 
84° 


90° 

89° 

90° 

90° 

90° 

90° 

90° 

90° 

80° 
81° 
82° 
84° 
86° 
87° 
89° 
90° 
90° 


78° 
80° 
81° 
83° 
85° 
86° 
88° 
89° 
90° 


77° 
79° 
80° 
82° 
84° 
86° 
87° 
89° 
90° 


For   Pg/y  =  c*  we  have  stability  for  all  0, 
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Uq/v    >  0.    For   0  <  Mq/u  1  t)(c  )  ,  max  b(c  )  ~  0.75, 


we  have  stability  for  0 


<  tt/2 


If  Mf^/M   increases,  a  belt  of  instability  appears, 
and  the  width  of  this  belt  increases  with  increasing  Mq/p 

In  the  table  below,  we  have  the  limits  of  the 
instability  belt   6,  and  Q^    —  unstable  for  Q-,  1  9  £  Q^  ' 

stable  otherwise,  for  some  values  of  c   and  y^/M* 

p      0 


Configuration  is  unstable  for  6,  1.  6  £  ^o 
(in  degrees) 


:^^- 

1. 

0 

1.25 

1.5 

2.0 

2.5 

3. 

0 

4.0 

5.0 

^1 

^2 

^1 

^2 

^1 

'2 

^1 

^2 

^1 

«2 

^1 

^2 

^1 

^2 

^1 

^2 

.2 

0 

0 

11 

28 

9 

37 

9 

47 

8 

52 

8 

56 

7 

60 

7 

63 

.3 

0 

0 

10 

29 

9 

39 

9 

49 

8 

54 

8 

58 

7 

62 

7 

65 

.4 

0 

0 

10 

32 

9 

41 

8 

51 

8 

57 

8 

60 

7 

64 

7 

67 

.5 

0 

0 

10 

34 

9 

44 

8 

59 

8 

59 

8 

63 

7 

67 

7 

69 

.6 

0 

0 

10 

37 

9 

46 

8 

57 

8 

62 

8 

66 

7 

70 

7 

72 

.7 

14 

17 

10 

39 

9 

49 

8 

59 

8 

65 

8 

68 

7 

72 

7  74 

.8 

13 

22 

10 

42 

9 

51 

8 

62 

8 

67 

' 

71 

7  75 

7  77 

.  .9 

12 

25 

10 

43 

9 

53 

8 

64 

8 

69 

1 

73 

7  77 

7  79 

1.0 

12 

28 

10 

44 

9 

54 

8 

65 

8  71 

1 

75 

7  79 

7  81 
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(c)   N  ^  4  . 

y./n  <  0 ,  for  small  values  of  | Vn/v \    we  have  a  polar 
cap  of  stability  which  decreases  with  increasing  values 


of  I  y^^/ij  I  .   For   yn/y  <  -1.0,  the  configuration  is  unstable 


in  the   range 


<  90° 


Given  below  are  the 


values  of  the  limit  6   of  the  cap  for  N  =  4,5,6,7, 


Values  of  e   for  N  =  4 
c 

* 

Configuration  is  stable,  9  £  9  <_  0 


<  5" 


■\yo/y 
"p\ 

-  0.1 

-0.25 

-  0.5 

-  1.0 

i 

.2 

48° 

37° 

10° 

<  5° 

.3 

47° 

36° 

9° 

<  5° 

.4 

45° 

33° 

9° 

<  5° 

.5 

44° 

31° 

9° 

<  5° 

.6 

42° 

28° 

8° 

<  5° 

.7 

39° 

26° 

8° 

<  5° 

.8 

37° 

24° 

8° 

<  5° 

.9 

35° 

22° 

8° 

<  5° 

1.0 

33° 

20° 

7° 

<  5° 

For  \ir./V'   =   c*    the  configuration  is  exponentially 
unstable. 


•124- 


Values  of  9   for  N  =  5 
c 

*  * 

Configuration  is  stable  for  e<e<e   ,6   <5' 


p   ^^ 

-  .1 

-0.25      -  0.5      -  1.0 

.2 

.3 
.4 
.5 
.6 
.7 
.8 
.9 
1.0 

39° 
38° 
37° 
35° 
34° 
32° 
31° 
29° 
28° 

31° 
30  = 
28° 
26° 
24° 
23° 
21° 
20° 
18° 

10° 
9° 
9° 
9° 
9° 
8° 
8° 
8° 
8° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

Values  of  0   for  N  =  6 
c 

* 

Configuration  is  stable  for  9  <_ 


-0.25 


-  0.5 


-  1.0 


.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 
1.0 


19° 
19° 
18° 
17° 
15° 
14° 
14° 
13° 
12° 


10° 

10° 

10° 

10° 

9° 

9° 

9° 

9° 

8° 


7° 
7° 
7° 
7° 
7° 
7° 
7° 
6° 
6° 


<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 
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Values  of  6   for  N  =  7 
c 

* 
Configuration  is  stable  for  0  <  e  < 


—   —  c 


<  5° 


4j^/m 

-  .1 

-0.25 

-  0.5 

-  1.0 

"P^-- 

.2 

10° 

8° 

7° 

<   5° 

.3 

10° 

8° 

7° 

<   5° 

.4 

9° 

8° 

7° 

<   5° 

.5 

90 

8° 

7° 

<   5° 

.6 

90 

8° 

7° 

<   5° 

.7 

90 

8° 

7° 

<   5° 

.8 

90 

8° 

7° 

<   5° 

.9 

90 

8° 

6° 

<   5° 

1.0 

90 

8° 

6° 

<   5°   ' 

Values  of  6   for  N  =  8 
c 

* 
Configuration  is  stable  for   6  <  6  < 


<    5° 


"\yo/y 

-  .1 

-0.25 

-  0.5 

-  1.0 

S^\ 

.2 

8° 

8° 

7° 

<  5° 

.3 

8° 

7° 

7° 

<  5° 

.4 

8° 

7° 

7° 

<  5° 

.5 

8° 

7° 

7° 

<  5° 

.6 

8° 

7° 

7° 

<  5° 

.7 

8° 

7° 

6° 

<  5° 

.8 

8° 

70 

6° 

<  5° 

.9 

8° 

7° 

6° 

<  5° 

1.0 

8° 

7° 

6° 

<  5° 

For   Prv/y   =    c*   we   have   exponential   instability 
for  N   =    4,5,6,7,8, 
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If  VIq/v  >  0 ,  we  have  a  cap  of  stability. 

As  Mf^/^l   increases,  this  cap  becomes  larger,  reaching 

a  point  where  we  have  stability  for  0  <_  9  £  Tr/2 .   Then  a 

belt  of  instability  appears.   With  increasing  Vq/ij,  the 

belt  of  instability  increases.   In  the  tables  given  below, 

we  have  the  limit  9   of  the  cap  and  the  limits  9w  Qj    °-^ 

the  belt  for  some  values  of  N,  y^/y,   c  . 

0      p 
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Values  of  6  . 
c 

Configuration  is  stable  for  0  _< 

N  =  4 


'"-^y 

p   ^ . 

.25 

.5 

1.0 

1.5 

.2 

66° 

75° 

90° 

90° 

.3 

66° 

76° 

90° 

90° 

.4 

67° 

77° 

90° 

90° 

.5 

67° 

79° 

90° 

90° 

.6 

69° 

81° 

90° 

90° 

.7 

70° 

82° 

90° 

90° 

.8 

72° 

84° 

90° 

90° 

.9 

74° 

86° 

90° 

90° 

1.0 

76° 

87° 

90° 

90° 

Values  of  Q  -.  ,    62 


Configuration  is  unstable  for  Q-i^O^G-;  stable  otherwise 


p--^o/^ 

c 

P 

- 

2.5 

3.0   1 

5 

0 

10.0 

11.0   1 

^1  '2 

^1 

'2 

h 

'2 

^1  '2 

^1 

'2 

.2 

13°  16° 

11° 

26° 

90 

42° 

8° 

51° 

8° 

52° 

.3 

13° 

15° 

11° 

26° 

90 

41° 

8° 

51° 

8° 

51° 

.4 

13° 

15° 

11° 

25° 

90 

40° 

8° 

50° 

8° 

51° 

.5 

- 

- 

11° 

24° 

90 

40° 

8° 

50° 

8° 

51° 

.6 

'  - 

11° 

24° 

90 

39° 

8° 

49° 

8° 

50° 

1   -7 

1  _ 

10° 

23° 

9° 

39° 

8° 

49° 

8° 

50° 

.8 

- 

10° 

22° 

90 

39° 

8° 

49° 

1  8° 

50° 

.9 

1  ~   1 

10° 

21° 

90 

38° 

8° 

49° 

7° 

50° 

1.0 

1  -   ■  -   1 

10° 

21° 

8° 

38° 

7° 

50°  '  '  7° 

50° 
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N  =  5 


Values  of 


Configuration  is  stable  for  0  _< 


<  e  , 

—  c 


^^q/*^ 

.25 

.5 

1.0 

1.5 

3.0 

s^--. 

1 

.2 

53° 

61° 

73° 

82° 

90°    i 

.3 

53° 

61° 

73° 

83° 

90° 

.4 

53° 

61° 

74° 

83° 

90° 

.5 

53° 

62° 

75° 

84° 

90°    ' 

.6 

52° 

62° 

76° 

84° 

90°    i 

.7 

52° 

63° 

76° 

84° 

90°    1 

.8 

52° 

64° 

77° 

85° 

90°    i 

i   .9 

53° 

64° 

78° 

85° 

90°     ; 

1.0 

53° 

65° 

78° 

85° 

90° 

Values  of  Q, ,    02  - 


Configuration  is  unstable  for  6,<_e<^02;  stable  otherwise, 


X^/i' 

- 

4.5 

5.0 

10.0 

11.0 

l=p  \ 

^1 

'2 

^1 

«2 

^1 

^2 

^1     '2 

1    -2 

!     .3 

1 

t     .4 

1     -^ 

!     .7 

.8 

.9 

1.0 

1 

13° 

15° 

11° 
11° 
11° 
12° 
12° 
12° 
13° 

21° 
20° 
18° 
17° 
16° 
14° 
13° 

_ 

9° 
9° 
9° 
9° 
9° 
9° 
9° 
9° 
8° 

38° 
37° 
35°  i 
34° 
32° 
30° 
29° 
27°  1 

26° 

1 

OOOOOOOOOOOOOOCOOO 
ooooooooo 

40° 
38° 

37° 
35° 
33° 
31° 
30° 
28° 
27°i 

1 

1 

1 
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N  =  6 


Values  of  6 


Configuration  is  stable  for  0  <_ 


<  e  . 

—  c 


.  p  _\> 

.2 

.25 

.5 

1.0 

1.5 

3.0 

5.0 

36° 

44° 

55° 

64° 

81° 

90° 

.3 

35° 

43° 

55° 

64° 

81° 

90°  i 

.4 

34° 

43° 

55° 

64° 

81° 

90° 

i   .5 

33° 

42° 

55° 

64° 

80° 

90° 

1 

.6 

32° 

41° 

54° 

64° 

80° 

90° 

!    .7 

31° 

40° 

54° 

64° 

79° 

90° 

.8 

30° 

40° 

54° 

64° 

79° 

88° 

J    ■' 

29° 

39° 

54° 

63° 

78° 

87° 

\     1.0 

28° 

38° 

54° 

63°   !   77° 

86° 

Values  of 


'1'  "2 


Configuration  is  unstable  for  ^-^±^±^2''    stable  otherwise. 


C  ^-x 

p 

7.5 

8    , 

10.0 

11.0 

^1 

'2 

'1 

'2 

^1 

'2 

^1   ^2 

.2 
.3 

.4 
.5 
.6 
.7 
.8 

•5   1 
1.0   1 

12° 
12° 
12° 

17° 
16° 

14° 

1 

11° 
11° 
11° 
11° 
11° 
\  12° 

1 

20° 
19° 
17° 
16° 
14° 

12° 

i 

10° 
10° 
10° 
10° 
10° 
10° 
10° 
10° 
10° 

27° 
25° 
23° 
22° 
20° 
18° 
17° 
15° 
14° 

9° 
90 

9° 
90 

90 

9° 
90 

90 

90 

29° 
27° 
25° 
23° 
22° 
20° 
18° 
17° 
16° 
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N 


Values  of  9 


Configuration  is  stable  for  0  £  9  <_  6 


>^ 

.25 

.5 

1.0 

1.5 

3.0 

5.0 

.2 

22° 

31° 

43° 

52° 

70° 

84° 

.3 

21° 

30° 

43° 

52° 

69° 

83° 

.4 

20° 

29° 

42° 

51° 

69° 

82° 

.5 

19° 

28° 

41° 

51° 

68° 

81° 

.6 

18° 

27° 

40° 

50° 

67° 

80° 

.7 

17° 

26° 

40° 

49° 

67° 

78° 

.8 

16° 

25° 

39° 

49° 

66° 

77° 

.9 

15° 

24° 

38° 

48° 

65° 

76° 

1.0 

14° 

23° 

37° 

47° 

64° 

74° 

Configuration    is   unstable    for    9, 


<_  Q    ;    stable  otherwise 


^p\ 

10. 

0 

11. 

0 

15. 

0 

\ 

'2 

^1 

«2 

^1 

'2 

.2 

12° 

13° 

10° 

17° 

9° 

25° 

.3 

- 

- 

10° 

16° 

9° 

24° 

.4 

- 

- 

11° 

14° 

9° 

22° 

.5 

- 

— 

11° 

13° 

9° 

20° 

.6 

- 

- 

11° 

12° 

9° 

18° 

.7 

- 

- 

- 

- 

9° 

17° 

.8 

- 

- 

- 

- 

9° 

15° 

.9 

89° 

90° 

- 

- 

9° 

14° 

1.0 

88° 

90° 

89° 

90°  ' 

9° 

13° 

-131- 


N  =  8 


Values  of 


Configuration  is  stable  for  0  < 


<  0  , 
—  c 


^--^/M 
S^\ 

.25 

.5 

1.0 

1.5 

3.0 

5.0 

10.0 

11.0 

.2 

10° 

13° 

27° 

38° 

57° 

72° 

90° 

90° 

.3 

10° 

13° 

26° 

37° 

56° 

71° 

90° 

90° 

.4 

10° 

12° 

25° 

36° 

55° 

70° 

88° 

90° 

.5 

10° 

12° 

24° 

35° 

54° 

68° 

86° 

88° 

.6 

10° 

12° 

23° 

34° 

53° 

67° 

84° 

86° 

.7 

10° 

12° 

22° 

33° 

52° 

66° 

82° 

84° 

.8 

9° 

11° 

21° 

31° 

51° 

64° 

80° 

81° 

.9 

9° 

11° 

20° 

30° 

50° 

63° 

78° 

79° 

1.0 

9° 

11° 

19° 

29° 

49° 

61° 

76° 

77° 

Configuration  is  unstable  for  6,  <^ 


<_   6_;  stable  otherwise. 


15 

20 

^1 

'2 

h 

'2 

'l 

'2 

^1 

^2 

.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 
1.0 

- 

10° 
10° 
10° 
10° 
10° 
10° 
10° 
86° 
84° 

16° 
15° 
14° 
14° 
13° 
12° 
11° 
90° 
90° 

9° 
9° 
9° 
9° 
9° 
9° 
9° 
9° 
9° 

23° 
22° 
20° 
18° 
17° 
15° 
14° 
13° 
13° 

- 
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A  Graphic   Representation   of   max    |Re    A. |    and   Im|A. |    for 
Different  Values    of   the    Parameters    are    Given 


Figure 
Number 

Plot 

Contents 
of 

5 

Versus 

For 

N= 

c 
P 

Mq/M 

A. 27 

max 

Re  X. 

1 

6 

3 

.2,  . . . , 

1.0 

1.0 

A. 28 

max 

Re  A. 

1 

e 

3 

.2,. . . , 

1.0 

1.5 

A. 29 

max 

Re  A. 

1 

e 

3 

.2,  .  .  .  , 

1.0 

2.0 

A. 30 

max 

Re  A. 

1 

e 

3 

.2,. ..  , 

1.0 

-4.0 

A. 31 

Im 

A. 

1 

e 

3 

1.0 

-4.0 

A. 32 

max 

Re  A. 

1 

0 

4 

.2, . .  ,  , 

1.0 

.5 

A. 33 

max 

Re  A^ 

e 

4 

.2,.  .  .  , 

1.0 

-  .5 

A. 34 

max 

Re  A. 

1 

e 

5 

.2,  .  .  .  , 

1.0 

.5 

A.  35 

max 

Re  A. 

e 

5 

.2,  .. .  , 

1.0 

1.0 

A. 36 

i 

max 

Re  A. 

e 

5 

.2,. . . , 

1.0 

1.5 

A. 37 

max 

Re  A. 

e 

6 

.2,. . . , 

1.0 

.5 
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A  Graphic  Representation  of  max  | Re  A. |  and  Im  |A. |  for 
Different  Values  of  the  Parameters  are  Given 

-  Continued  - 


Figure 
Number 

Contents 
Plot  of 

Versus 

For 

N  = 

c 

P 

Pq/M 

A. 38 

max   Re  A 

0 

6 

.2,. 

..,  1.0 

1.0 

A.  39 

max   Re  A . 

1 

6 

6 

.2,. 

..,  1.0 

1.5 

A. 40 

max   Re  A . 

1 

6 

6 

.2,  . 

.. ,  1.0 

-  .5 

A. 41 

max   Re  A . 

e 

7 

.2,. 

.. ,  1.0 

.5 

A. 42 

max   Re  A . 

e 

7 

.2,  . 

. .,  1.0 

1.0 

A. 43 

max   Re  A . 

1 

e 

7 

.2,  . 

. . ,  1.0 

-  .5 

A.  44 

max   Re  A . 

1 

6 

3 

.2,  . 

.. ,  1.0 

.5 

A. 45 

max   Re  A . 

e 

8 

.2,. 

.  .  ,  1.0 

1.0 

A. 46 

max   Re  A . 

1 

9 

8 

.2,  . 

. . ,  1.0 

1.5 

A. 47 

max   Re  A . 

0 

8 

.2,  . 

. .,  1.0 

-  .5 
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Huo/mj:  1 .0000 


.CP.  li) 


_^ct^^2A£-^ 


°i.m  II  00      2100     79.00      n.oo      ^S.OO      SI-OO      bl.OO 

COLHIITUnC 


MUO/MUr    1  .SOCO 


tC.slJ) 


•^,00  p. 00  ji.oo  <i.oo  5). 50  w.oo  n.ofl  BI.OO 

COLHIITUDE 


Figure  A. 27 


Figure  A. 2 8 


Ht 


n=3 


Mix)/MU=  e.oooo 


^.00  n.oQ         w.oo         ti.oo         53.00         «.oa         ti.oo         of.on 

COLflTlTUDE 


Figure   A. 29 


l«X>/MU=-<I.OOOO 


c.r. 


i 

■7 
■8 


IT.OO  O.OO  tl.W  n.OO  8».0O  TT.CU  w.oo 

COLflTirUDE 


Figure   A. 30 
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C.P.ct.O 
t«Xl/Mlte-M  .0000 


17.00  tf.OO  41. OO  O.OO  fH.n  77.00  69.00 

COLfillTUDE 


MUO/HU=      .5000 


S 


8  


».aO  It.OO  53.00 

COLflTiruDE 


6E.ra  7V.0U  f;9.00 


Figure  A. 31 


Figure   A. 32 


HUO/MUr    -.5000 


77.00  Bt.OO 


8 

- 

CI! 

I 

f- 

N=5 

/,.. 

S 

( 

MUO/MU= 

.SOOO 

/ 

» 

/-•^ 

o 

/ 

.5 

or. 

.6 

i« 

y^ 

1 

.^ 

t 

^  .8 

s 

^.-^  -9 

i 

i     /X 

8 

i 

"^.OO  17.00 


Figure   A. 33 


79. CO  <41.0D  SJ.OO 

MLfllllUDE 

Figure   A. 34 


6S.0O  77.00  e».M 
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N=S 


MUO/  MU.l  .0000 


*^.00  20.00 


40.00  63.00  ao.oo 

COLfllllUDE 

Figure   A. 35 


lOO.OO  )».00  HO  .00 


MUO/MU=    1 .5000 


1 


.e.jLO 


17.00  M.eo  III. 00  SI. 00 

counniuoE 


u.oo  IT. 00  el.oc 


Figure   A. 36 


B 

: 

N=5 

•0 

MUO/MU= 

.5000 

cx 

1 

.2 

5 

/C— -■• 

„ 

//■ 

T 

o 

/C—-^ 

^'  ■ 

!/^ 

/^                                      ■ -6 

s 

r^^^^^-T^^^^^"  ■' 

8 

ffli           ^^^-^         ~^ — —  'S 

S 

j 

^--^^.0 

8 

1 

F 

,              / 

1 
I. 

^.00  17.00  ?9.0O  41.00  63.00 

COlflTIIllOE 

Figure  A. 3  7 


66.00  77.00  61.00 


8 


MUO/MU;    1 .0000 


6.00  17.00 


Oi.OO  41.00  63  GU 

COLRIiruOE 

Figure   A. 38 


es.oo  77 .00  St. 00 


-137- 


Nr6 

nuo/hu= 

1  .sooo 

C.f. 

i 

/v^l 

W^"' 

f^ 

8 


rwo/iiUr  -.5000 


2t  N  41. n  S3  00  U 

COLHI 1 TUOt 

Figure   A. 39 


?g.00  Mt.OO  53.00  6! 

COLBTITUOE 

Figure   A. 40 


%M 


Figure   A. 41 


29.00  m.OO  53.00 

COLBIltUDE 

Figure   A. 42 
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^.ai 


»1U0/hU=  -  .SflOO 


29.00  41.00  SJ.OD 

COlfilltUOt 


Figure   A. 43 


66.00  11.00  99.00 


«uo/nu=     .5000 


*^.0«  17.00 


21. OO  11.00  S9.00 

roiniiTuoE 


Fiaure   A. 44 


77.00  •9.00  . 


"^.00 


11.00  29.00  HI  .00  51.00 

COLflTlTUOE 


65. 00      71.00      69.00 


*^.00       n.OO      M  OU      HI. 00 53-00 

LOLRTITUOE 


ftS.OO  77.00  «■» 


Figure   A. 45 
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Figure   A. 46 


MUO/KU-    -.5000 


29.00  11. OU  ^i.OU 

COlflll'UOt 


6^.00  77-00 


Figure   A. 4  7 
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(iv)   N  Vortices  on  a  Circle  of  Colatitude,  One  Polar  Vortex, 
Boundary  Condition   at   8  =  Tr/2 .    (2  <  N  <  8) 


(a) 


N  -  2 


Mq/v    >    0, 


We  have  a  belt  of  stability,  the  limits  6,  and  e„  of  the 


belt  are  given  below  for  some  values  of  c    and  Un/u 

p       0 


Configuration  is  stable  for  6^  <_ 


'2- 


(in  degrees) 

0.125 

.25  1 

.5 

1.0 

1.5 

3.0 

4.0 

^1 

'2 

^1 

'2 

^1 

^2 

^1 

'2 

^1  '2 

^1 

'2 

^1 

'2 

.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 
1.0 

65 
65 
65 
65 
65 
65 
65 
64 
64 

66 
66 
66 
67 
67 
68 
69 
69 
70 

67 
67 
67 
67 
67 
66 
66 
66 
66 

70 
71 
71 
71 
72 
72 
73 
73 
74 

70 
69 
69 
69 
69 
68 
68 
68 
67 

76 
76 
70 
77 
77 
77 
77 
77 
78 

72 
71 
71 
71 
70 
70 
70 
69 
69 

82 
82 
82 
82 
81 
81 
81 
81 
81 

72 
72 
72 
72 
71 
71 
70 
70 
70 

84 
84 
84 
84 
84 
83 
83 
83 
83 

71 
71 
72 
72 
72 
72 
71 
71 
71 

>85 
>85 
>85 
>85 
>85 
>85 
>85 
>85 
>85 

63 
63 
70 
71 
71 
71 
71 
71 
71 

>85 
>85 
>85 
>85 
>85 
>85 
>85 
>85 
>85 

For    3.5    <_  Mrv/y    <    4.0      another   belt   appears   between 


60°  and  70°  colatitude.   The  values  for  this  belt  are  not 
tabulated.   It  can  be  seen  in  Figure  A. 69  a. 
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Pq/p  <  0. 

For  -1.0<_  iJp./iJ  <  0   we  have  instability.  For  Vq/m  -    -1.25 
a  belt  of  stability  appears  and  for  -2  <_  y^j/u  £  -1.5   we  have 
a  cap  of  stability. 

The  limits  9^  and  0   of  the  belt  of  stability  and 
the  limit  0   of  the  cap  are  given  below  for  various  values 
of  c   and  Pq/m. 


Configuration  is  stable  for  6,  <_ 


)^;    and  e  <_ 


■<5' 


P    \ 

-1.25 

-1.5 

-2.0 

^1 

'2 

e 

c 

9 
c 

.2 

11° 

47° 

57° 

65° 

.3 

11° 

45° 

55° 

63° 

.4 

12° 

43° 

53° 

61° 

.5 

13° 

42° 

51° 

58° 

.6 

18° 

40° 

50° 

57° 

.7 

26° 

39° 

49° 

55° 

.8 

32° 

37° 

48° 

54° 

.9 

unstable 

48° 

54° 

1.0 

unstable 
1 

48° 

53° 

For  Mf^/y  <  -2.5  ,  this  cap  is  broken  into  a  smaller 
cap  and  a  belt  of  stability.  In  the  tables  below  we  have 
the  limits  of  the  cap  and  belt  for  different  values  of  c 

IT 

and  Vq/u    ' 
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^q/u  =  -2.5 
Configuration  is  stable  for  6^  1  9  1  9^ 


c 

p 

^1 

'2 

^1 

^2 

.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 
1.0 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

35° 
31° 
26° 
63° 
61° 
59° 
58° 
57° 
56° 

45° 
37° 
28° 

69° 
67° 
65° 

1 

Mq/m  =  -3.0 


Configuration  is  stable  for  6,  <_ 


<  ^2 


c 

p 

^1 

62 

^1 

'2 

.2 

<  5° 

47° 

60° 

72° 

.3 

<  5° 

43° 

55° 

70 

.4 

<  5° 

39° 

47° 

68° 

.5 

<  5° 

34° 

38° 

66° 

.6 

<  5° 

64° 

- 

- 

.7 

<  5° 

62° 

- 

- 

.8 

<  5° 

61° 

- 

- 

.9 

<  5° 

60° 

- 

- 

1.0 

<  5° 

59° 

- 

- 
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Pq/m  =  -3.5 
Configuration  is  stable  for  9,  £  6  S  ^o 


c 

p 

^1 

^2      ^1 

^2 

.2 

<  5° 

54° 

67° 

74° 

.3 

<  5° 

50° 

64° 

72° 

.4 

<  5° 

46° 

59° 

70° 

.5 

<  5° 

42° 

51° 

68° 

.6 

<  5° 

37° 

42° 

66° 

.7 

<  5° 

33° 

35° 

64° 

.8 

<  5° 

63° 

- 

- 

.9 

<  5° 

61° 

- 

- 

;  1.0 

<  5° 

61° 

- 

- 

Vq/V    =    -4.0 
Configuration  is  stable  for  6,  1  6  1.  Q^ 


"p 

^1 

'2 

^1 

'2 

.2 

<  5° 

59° 

72° 

76° 

.3 

<  5° 

55° 

69° 

74° 

.4 

<  5° 

51° 

65° 

72° 

.5 

<  5° 

47° 

60° 

70° 

.6 

<  5° 

43° 

53° 

68° 

.7 

<  5° 

39° 

45° 

66° 

.8 

<  5° 

35° 

38° 

64° 

.9 

<  5° 

63° 

- 

- 

1.0 

<  5° 

62° 

- 

- 
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Uq/M  =  -4.5 
Configuration  is  stable  for  9,  <_  Q    1  Qj 


c 

p 

'l 

«2 

\ 

^2 

.2 

<    5° 

63° 

74° 

77° 

.3 

<  5° 

59° 

72° 

76° 

.4 

<  5° 

55° 

69° 

74° 

.5 

<  5° 

51° 

65° 

71°  I 

.6 

<  5° 

48° 

60° 

69° 

.7 

<  5° 

44° 

54° 

67°  , 

.8 

<  5° 

41° 

47° 

66° 

.9 

<  5° 

37° 

40° 

64° 

1.0 

<  5° 

33° 

35° 

63° 

Mq/M  -  -5.0 
Configuration  is  stable  for  6,  £  6  S  Q-, 


c 

p 

^1 

«2 

'i 

«2 

.2 

<  5° 

65° 

76° 

78° 

.3 

<  5° 

62° 

74° 

77° 

.4 

<  5° 

58° 

71° 

75° 

.5 

<  5° 

55° 

68° 

73° 

.6 

<  5° 

51° 

64° 

70° 

.7 

<  5° 

48° 

60° 

68° 

.8 

<  5° 

45° 

54° 

67° 

1 

.9 

<  5° 

42° 

47° 

65° 

1.0 

<  5° 

38°  1   41° 
1 

64° 
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\1q/v    =    -6.0 
configuration   is    stable    for    9,    <_  0    1  02 


"p 

^1 

^2 

^1 

'2 

.2 

<  5° 

70° 

79° 

— t 

80° 

.3 

<  5° 

67° 

77° 

78° 

.4 

<  5° 

63° 

75° 

77° 

.5 

<  5° 

60° 

72° 

75° 

.6 

<  5° 

56° 

69° 

72° 

.7 

<  5° 

53° 

66° 

70° 

.8 

<  5° 

50° 

63° 

69° 

.9 

<  5° 

48° 

59° 

67° 

1.0 
1 

<  5° 

46° 

54° 

66° 
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N    =    3 


(a)      \1q/]i    <   0. 

In    this   case   we   have    a   cap   of   stability. 

Given   below   are    some   values   of    9    . 

c 

Value  of  6  .  Configuration  is  stable  for  0  <_ 


<    6  , 
—  c 


^o/p 

-0.25 

-  .5 

-1.0 

-1.5 

-2.0 

-2.5 

s^\ 

.2 

43° 

42° 

39° 

35° 

30° 

23° 

.3 

43° 

42° 

39° 

35° 

30° 

23° 

.4 

43° 

42° 

39° 

35° 

31° 

24° 

.5 

43° 

42° 

39° 

36° 

31° 

24°   1 

.6 

44° 

42° 

40° 

36° 

32° 

24° 

.7 

44° 

43° 

40° 

37° 

32° 

25° 

.8 

44° 

43° 

40° 

37° 

32° 

25° 

.9 

44° 

43° 

41° 

37° 

33° 

26° 

1.0 

44° 

43° 

41° 

38° 

34° 

26° 

Value  of 


Configuration  is  stable  for  0  < 


c 

p 

\lQ/Vi  =  C* 

-3.0 

-4.0 

-6.0 

.2 

39° 

45° 

63° 

72° 

.3 

39° 

41° 

59° 

70° 

.4 

40° 

37° 

55° 

66° 

.5 

41° 

35° 

52° 

63° 

.6 

42° 

34° 

50° 

60° 

.7 

43° 

34° 

48° 

57° 

.8 

43° 

34° 

47° 

56° 

.9 

44° 

35° 

46° 

54° 

1.0 

44° 

36° 

46° 

53° 
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(b)   y^/y  >  0. 

For   0<  Vq/v    <  .85,  we  have  a  cap  of  stability  defined 

by  0  <  6   .   In  the  table  below  we  have  the  values  of  0 
-^    —  c  c 

for  Uq/p  =  -5,  .75,  .85. 


Values  of 


Configuration  is  stable  for  0  <  0  < 


0.25 

.5 

.75 

.85 

.2 

46° 

47° 

48° 

49° 

.3 

46° 

47° 

48° 

48° 

.4 

46° 

47° 

48° 

48° 

.5 

46° 

47° 

48° 

48° 

.6 

46° 

47° 

48° 

48° 

.7 

46° 

47° 

47° 

48° 

.8 

46° 

47° 

47° 

47° 

.9 

46° 

47° 

47° 

47° 

1.0 

46° 

47° 

47° 

47° 

For  larger  values  of  V^q/v^,       -85  <  Uq/ij  <  1-5,  this 
cap  of  stability  is  broken  into  a  smaller  cap  and  a  belt 
of  stability. 

In  the  table  below  we  have  the  limits  of  the  smaller 
cap  and  the  belt  of  stability  for  values  of  Mq/ij  =  1.0  and  1.5 
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Configuration  is  stable  for  0  <_  G  £  9   and 


I,  <_  e  <_  62. 


^P " 

Mq/M  =  1- 

0 

Uq/u   =   1. 

5 

cap  lim 

limits 

of  belt 

cap  lim 

limits 

of  belt 

^c 

^1   ■ 

■    '2 

'c 

^1  1 

-   '2 

.2 

49° 

- 

- 

8° 

43° 

52° 

.3 

49° 

- 

- 

8° 

45° 

51° 

.4 

49° 

- 

8° 

47° 

50° 

.5 

49° 

- 

- 

8° 

- 

.6 

48° 

- 

- 

8° 

- 

- 

.7 

13° 

18° 

48° 

8° 

- 

- 

.8 

12° 

23° 

48° 

8° 

- 

- 

.9 

11° 

27° 

47° 

8° 

- 

- 

1.0 

11° 

29° 

47° 

8° 

- 

- 

The  belt  disappears  for  increasing  values  of  Mq/m 
and  for  M^x/y  =  2.0,  2.5;  we  have  again  a  cap  of  stability. 

Values  of  6  .   Configuration  is  stable  for  0  <_  6  <_  9  . 
c  c 


c 

p 



Pq/P  =  2.0 

Mg/y  =  2.5 

.2 

8° 

7° 

.3 

8° 

7° 

.4 

7° 

70 

.5 

70 

70 

.6 

7° 

7° 

.7 

7° 

7° 

.8 

70 

70 

.9 

70 

70 

1.0 

70 

70 
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Increasing  Mq/M/  the  new  cap,  is  broken  into  a  smaller 
cap  and  a  belt  of  stability.   Given  below  are  some  values  of 
the  limit   6   of  the  cap  and  Q-i/Qo  °^   ^^®  belt   for 
different  values  of  Mf^/^l• 


Configuration  is  stable  for   0  <_  0  <_  6 


'l  1  9  1  92- 


yg/y 

3.0 

4.0 

6.0 

c 

p 

cap 

be] 
^1 

Lt 
'2 

cap 

e 

c 

be] 

^1 

Lt 
^2 

cap 

6 
c 

be: 

^1 

Lt 
'2 

.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 
1.0 

7° 
7° 
7° 

7° 

1 

1  7° 

7° 

7° 

i 
1 
I  7° 

1 

7° 

60° 
63° 
64° 

63° 
63° 
64° 

6° 
6° 
6° 
6° 
6° 
6° 
6° 
6° 
6° 

63° 
63° 
62° 
65° 
66° 
67° 
66° 
67° 
68° 

67° 
67° 
68° 
68° 
68° 
68° 
67° 
68° 
68° 

6° 
6° 
6° 
6° 
6° 
6° 
6° 
6° 
6° 

67° 
67° 
66° 
66° 
65° 
65° 
67° 
68° 
68° 

73° 
73° 
73° 
73° 
72° 
72° 
72° 
72° 
71° 
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N  >  4 


Mq/P  <  0. 

A  cap  of  stability  is  again  common  for  n  >^  4 .   The  limits 


of  the  cap,  G   ,  are  given  below  for  N  =  4,5,6,7,8  and 


C       •  ^  f   •  -^  / 


.,  1.0  for  different  values  of  Vq/v • 


N   =    4 


Values    of    6    .      Configuration    is    stable    for   0    <    6    <    6    . 
c  —        —     c 


N=5 


p   \ 

-  .1 

-0.25 

-  0.5 

-  1.0 

.2 

.3 
.4 
.5 
.6 
.7 
.8 
.9 
1.0 

39° 
39° 
38° 
38° 
37° 
36° 
35° 
34° 
32° 

33° 
32° 
31° 
29° 
28° 
26° 
24° 
22° 
20° 

10° 
9° 
9° 
9° 
8° 
8° 
8° 
8° 
7" 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

<  5° 

Values  of 


Configuration  is  stable  for  0  <_ 


-  .1 

-0.25 

-  0.5 

-  1.0  1 

.2 

35° 

29° 

10° 

5° 

.3 

35° 

28° 

9° 

5°  ! 

1 

.4 

34° 

27° 

9° 

5°  ' 

.5 

33° 

26° 

99 

5° 

.6 

33° 

24° 

90 

5° 

.7 

31° 

23° 

8° 

<  5° 

.8 

30° 

21° 

8° 

<  5° 

.9 

29° 

20° 

8° 

<  5° 

1.0 

28° 

18° 

8° 

<  5° 
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N  =  6 


Values  of  6  .   Conf iqurationis  stable  for  0  < 
c  — 


\yo/M 

-  .1 

-0.25 

-  0.5 

-  1.0 

p   \ 

.2 

19° 

10° 

7° 

5° 

.3 

18° 

10° 

7° 

5° 

.4 

17° 

10° 

7° 

5° 

■= 

16° 

10° 

7° 

5° 

.6 

15° 

go 

7° 

5° 

.7 

14° 

9° 

7° 

5° 

.8 

14° 

go 

7° 

5° 

.9 

13° 

9° 

6° 

5° 

1.0 

12° 

8° 

6° 

5° 

N    =    7.      Values    of 


Configuration   is    stable    for   0    < 


<    6     , 
—      c 


^\J^o/^ 

P   ^^ 

-  .1 

-0.25 

-  0.5 

-  1.0 

.2 

10° 

8° 

7° 

5° 

.3 

10° 

8° 

7° 

5° 

.4 

go 

8° 

7° 

5° 

.5 

go 

8° 

7° 

5° 

.6 

go 

8° 

7° 

5° 

.7 

go 

8° 

7° 

5° 

.8 

go 

8° 

7° 

5° 

.9 

go 

8° 

6° 

5° 

1.0 

go 

8° 

6° 

5° 
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N  =  8.   Values  of  6  .   Configuration  is  stable   for  0  <_  6  <_  6^ 


P    ^^ 

-  .1 

-0.25  , 

- 

1 
-  1.0 

.2 

8° 

8° 

7°   ; 

j 

5° 

.3 

8° 

7° 

7° 

5° 

.4 

8° 

7° 

7° 

5° 

.5 

8° 

7° 

7° 

5° 

.6 

8° 

7° 

7° 

5° 

.7 

8° 

7° 

6° 

5° 

.8 

8° 

7° 

6° 

5° 

.9 

8° 

7° 

6° 

5° 

1.0 

8° 

7° 

6° 

5° 
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^q/U  >  0. 


Here  we  have  a  cap  of  linear  stability,  as  y^/y 
increases,  this  cap  is  broken  into  a  smaller  cap  and  a  belt 
of  stability.   Below  are  given  the  limits  9   of  the  cap, 


the  limit  9-,  of  the  smaller  cap  and  the  limits  82'^ 2 

nt  values  of  N,  u^/y  and  c  . 

0        p 

Configuration  is  stable  for  0  <_  9  £ 


of  the  belt  for  different  values  of  N,  Mq/m  and  c 


N  =  4.   Values  of 


""^q/^ 
S^\ 

.25 

.5 

1.0 

1.5 

1 

1 

.2 

47° 

50° 

52° 

54° 

.3 

470 

50° 

52° 

54° 

.4 

47° 

50° 

52° 

54° 

.5 

48° 

50° 

52° 

54° 

.6 

48° 

50° 

52° 

54° 

.7 

48° 

50° 

52° 

54° 

.8 

49° 

51° 

52° 

54° 

.9 

49° 

51° 

52° 

54° 

1.0 

50° 

51° 

52° 

54° 

Configuration  is  stable  for 


0    ±    Q    ±   Qy 


)<^9o;  unstable 

otherwise. 


t,' 

3.0       1 

5.0 

^1 

'2 

'3 

^1 

'2 

S 

.2 

10° 

27° 

59° 

12° 

17° 

58° 

.3 

10° 

26° 

59° 

12° 

16° 

58° 

.4 

10° 

26° 

59° 

12° 

15° 

57° 

.5 

10° 

25° 

59° 

57° 

- 

- 

.6 

10° 

24° 

58° 

57° 

- 

- 

.7 

90 

24° 

58° 

57° 

- 

- 

.8 

90 

23° 

58° 

57° 

- 

- 

.9 

9° 

22° 

58° 

56° 

- 

- 

1.0 

90 

22° 

57° 

56° 

- 

- 
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Configuration   is    stable    for   0    <_  Q    <_  Q,,    00-^-^3    ' 

unstable  otherwise. 


\^q/V 

5.0 

10.0 

11.0 

c 

p 

cap 

be] 

'2 

Lt 
^3 

cap 

^1 

be] 

'2 

Lt 
^3 

cap 

^1 

belt 

02    63 

.2 
.3 
.4 
.5 
.6 
.7 
.8 
.9 
1.0 

8° 
8° 
8° 
8° 
8° 
8° 
8° 
8° 
7° 

1 
41° 

41° 

41° 

40° 

40° 

40° 

40° 

39° 

39° 

64° 
64° 
63° 
63° 
63° 
62° 
62° 
61° 
61° 

7° 
7° 
7° 
7° 
7° 
7° 
7° 
7° 
6° 

51° 
51° 
51° 
51° 
51° 
51° 
51° 
51° 
51° 

70° 
69° 
69° 
68° 
68° 
67° 
67° 
66° 
66° 

7° 
7° 
7° 
7° 
7° 
7° 
7° 
6° 
6° 

52° 
52° 
52° 
52° 
52° 
52° 
52° 
52° 
52° 

70° 
70° 
70° 
69° 
69° 
68° 
67° 
67° 
66° 

N  =  5 


Values  of 


Configuration  is  stable  for  0  <_ 


<  e  , 
—  c 


p  ^ 

.25 

.5 

1.0 

1.5 

3.0 

.2 

44° 

48° 

53° 

55° 

59° 

.3 

44° 

48° 

53° 

55° 

59° 

.4 

44° 

48° 

53° 

55° 

59° 

.5 

44° 

49° 

53° 

55° 

59° 

.6 

45° 

49° 

54° 

55° 

59° 

.7 

45° 

49° 

54° 

55° 

59° 

.8 

45° 

50° 

54° 

55° 

59° 

.9 

46° 

50° 

54° 

55° 

59° 

1.0 

46° 

51° 

54° 

56° 

58° 
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Configuration    is    stable    for   0    <_  9    <_  6    ,    9;>   £  6   £   6^    ; 

unstable   otherwise. 


yg/iJ 

4.5 

5.0 

10.0 

11.0 

-- 

c 

p 

cap 

belt 

cap 

belt 

cap.  belt 

cap 

belt 

^1 

^2   ^3 

Ol 

^2 

1 
^3 

6,  e^ 

S 

^1 

^2 

^3 

.2 

.3 

.4 

.5 

.6 

.7  i 

.8 

1.0 

12° 
62° 
62° 
62° 
62° 
61° 
61° 
61° 
61° 

16° 

63° 

10° 
10° 
10° 

11° 

11° 
11° 

12° 
62° 
61° 

21° 
20° 
19° 
18° 
17° 
15° 
14° 

64°    8° 
63°    8° 
63°    8° 
63°    8° 
62°    8° 
62°    8° 
62°    8° 
8° 
7° 

38° 
37° 
36° 
34° 
33° 
31° 
29° 
28° 
27° 

69° 
69° 
69° 
68° 
68° 
67° 
67° 
66° 
66° 

7° 
7° 

7° 

1 
7° 

7° 

i  7° 

;   7° 

1 

.  7° 
7° 

39° 
38° 
37° 
35° 
34° 
32° 
31° 
29° 
28° 

70° 
70° 
69° 
69° 
68° 
68° 
67° 
67° 
66° 

N  =  6 


Values  of 


Configuration  is  stable  for  0  < 


<  e  , 

—  c 


1. 


.25 

34° 
33° 
33° 
32° 
31° 
31° 
30° 
29° 
28° 


.5 

40° 
40° 
39° 
39° 
39° 
39° 
38° 
38° 
38° 


1.0 

47° 
48° 
48° 
48° 
48° 
48° 
49° 
49° 
49° 


1.5 

52° 
52° 
53° 
53° 
53° 
53° 
53° 
54° 
54° 


3.0 

59° 
59° 
59° 
59° 
59° 
59° 
59° 
59° 
59° 


5.0 

63° 
63° 
63° 
63° 
62° 
62° 
62° 
62° 
62° 
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Configuration  is  Stable  for  0  <  9  <  G,;  Q^    <_ 

Unstable  otherwise. 


Pq/^ 

7.5 

10.0 

1 

11.0 

c 

p 

cap 

belt 

cap 

belt 

cap 

belt 

^1 

^2 

^3 

^1 

^2   ^3 

^1 

'2      '3 

.2 

11° 

18° 

67° 

9° 

27° 

69° 

8° 
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Configuration   is    Stable    for   0    <_  G    £  Q-i 
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Configuration  is  stable  for  0  _<  6  <_  9   ,  ^t  -  ®  -  ^3 
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is   given  for  different  values  of  the  parameters. 
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By  using  the  tables,  we  can  construct  the  curves  of 
linear  stability  which  separate  the  unstable  from  the 
stable  regions,  in  the  plane  (y  /y,  0)   for  each  N  and  c  . 
In  Figures  A.69a,b,c  and  A. 70,  the  curves  are  drawn  for 
N=2,...,8,   c  =.2   with  and  without  boundary  conditions 
respectively . 

Comparing  the  figures,  we  can  see  that  the  effect  of 
the  boundary  condition   is  to  reduce  the  zones  of  stability. 
Another  feature  is  that  for  smaller  colatitudes,  N  >^  4 , 
G  _<  40°,   Mrj/y  >  3,   the  zones  of  instability   are  almost 
the  same  with  and  without  boundary  conditions,  which  means 
that  in  the  region  mentioned,  the  stability  properties  do 
not  depend  on  the  existence  of  the  boundary  conditions. 
Also  for  small  colatitude ,  the  coriolis  parameter  does  not 
play  a  significant  role  in  the  stability  properties.  Both 
phenomena  were  present  when  considering  the  cases  with 
no  polar  vortex. 

Note  that  for  N  =  2,  with  and  without  boundary  conditions, 
lif./v   =  0   is  in  the  zone  of  instability. 

For  the  cases  with  polar  vortex,  two  families  of  eigen- 
values are  evident,  one  family  is  pertinent   to  the  circle 
vortices  and  yield   the  left  branches  of  the  stability 
curves,  the  other  family  is  pertinent  to  the  polar  vortex 
and  yields  the  right  branches  of  the  stability  curves. 

For  the  case  of  no  pole,  and  no  boundary  condition,   we 
had   f2  =  0  at  9  =  Tr/2 ,  and  this  introduced  another  pair  of 
zero  eigenvalues. 
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When  a  pole  is  added,  0,^0      at   6  =  tt/2,  and  the 
singularity   is  removed  from   6  -   Tr/2   by  the  effect  of  the 

polar  vortex  and  placed  at  the  value  of  9  such  that 

* 

]X^/W    =    C  (6,N,C  )  . 

When  there  is  no  polar  vortex,  the  effect  of  the 
coriolis  parameter  can  be  summarized  as  follows: 

(a)  No  boundary  condition: 

For  a  given  configuration,  increasing  the  coriolis 
parameter   c   ,   the  polar  cap  of  stability  decreases 
for  all  N. 

(b)  Boundary  Condition 

For  a  given  configuration,  increasing   the  coriolis 
parameter,  the  polar  cap  of  stability  increases  for  N  =  2,3 
and  decreases  for  N  >  3. 

When  a  pole  is  present  there  is  no  such  simple  rule. 

For  example,  for  N  =  2,  no  boundary  condition, 

Vr./V'    <  0  ,  we  have  a  belt  of  stability,  if  we  increase  c   , 

the  belt  decreases.   For  the  same  case  but  1Jq/1-i  >  0,  the 

belt  of  instability  increases,  i.e.   stability  decreases, 

when    c    increases.   For  N  -  4,  no  boundary  condition, 
P 

U-/y  <  0,  the  polar  cap  of  stability  decreases  with  increasing 
c  .   For  N  =  5,  boundary  condition,  V-q/v   =  -25,  the  cap  of 

It 

stability   increases   with   increasing      c      but    for  N   =   6, 
boundary   condition,     Vrv/P  =    -25,    the   cap   of   stability   decreases 

with   increasing    c    • 

P 
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2 .    Numerical  Solutions  of  the  Nonlinear  Equations  of 
Motion  for  Configurations  with  No  Polar  Vortex, 
Both  with  and  Without  Boundary  Condition. 

In  Section  3,  we  have  discussed  configurations  with 
polar  vortex  and  given  a  description  of  the  motions  in  the 
zones  of  exponential  stability  and  instability. 

For  configurations  without  polar  vortex,  we  have  seen 
that  if   c   =  0,  there  are  no  linearly  increasing  solutions; 
whereas   c,  /  0  implies  that  the  longitudinal  components 
increase   linearly  with  time. 

Configurations  consisting  of  four  vortices  on  a  circle 
of  colatitude  are  studied  in  this  section.   N  =  4  was 
chosen  since  it  is  the  smallest  number  of  circle  vortices 
with  configurations  in  the  exponentially  stable  and  unstable 
zones   for  both  cases,  with  and  without  boundary  conditions. 

In  the  figures  to  be  given  later  on,  we  have  not  plotted 
(J),  (t),  k  =  1,...,4   but  instead 

D^(t)  =  4)^  -  (*i(t)  +  (N-1)  ^  ] 

Dj^(t)  =  (^j^(t)  -  (4)j^_i(t)  +   ^}     '         k  =  2,3,4. 

D  (t)  gives  a  measure  of  the  relative  oscillations, 
in  the   (J)   direction,  between  the  (k-l)st  and  k-th  vortex. 
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(i)  4  Vortices  with  Equal  Strength  y ,  Equally  Spaced  on 
a  Circle  of  Colatitude  9 ,  No  Boundary  Condition  at 
the  Equator. 

Here  we  study  three  cases  for  the  coriolis  parameter 

c   =  .7   and  one  case  for  the  coriolis  parameter  c  =  .5. 
P  P 

The  linear  analysis  establishes  that  zero  is  an  eigenvalue 
of  multiplicity  two  for  all  values  of  the  parameters. 

In  particular   for  c  =    .1 ,      we  find  that  for 
6  <_  49.6864°,  the  configuration  is  exponentially  stable; 
whereas  it  is  exponentially  unstable  for  6  >_  49.6865°. 

For  c  =  .5,  the  configuration   studied  is  in  the 
exponentially  ions  table  region. 

(a)   9  =  45°,  c  =  .7.   The  configuration  is  exponentially 
stable. 

The  experiments  were  carried  out  with  initial  displace- 
ments  eg,,  ea,   between  +1°  and  +  30°.   The  results  show 
that  all  the  resulting  motions  are  linear  stable  motions, 
even  though   in  some  cases,  the  amplitudes  of  the  resulting 
oscillations  are  bounded  by  3/2  of  the  initial  disturbance. 
This  is    particularly  true  when  the  initial  displacement 
is  in  the  longitudinal  direction. 

For  initial  displacements  in  the   G   direction,  the 
circle  vortices  tend  to  oscillate  about  a  new  position  north 
or  south  of  the  equilibrium  circle  of  colatitude   depending 
on  the  sign  of  eS-,-   This  phenomenon  has  the  effect  of 
changing  the  angular  velocity  of  rotation. 
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In  the  table  given  below,  we  have  listed  nonlinear 
periods  of  rotation,  for  different  types ^^d  magnitudes  of 
the  initial  displacement. 


Configurations  at  45°  Colatitude,   c   =  .7,  \i   =    .15*  10 

p 

Equilibrium  Period:   677.66  [hours] 

Nonlinear  Period  Initial  Displacement 
(average) 

707.344  [hours]  eB^  =  1° 

722.894  eB-,  =  2° 

737.723  eBj^  =  3° 

770.285  eB^  =  5° 

648.181  eB-,  =  -3° 

636.607  eB,  =  -5.0° 

677.607  ea   =   1° 

669.103  ea  =   5° 

682.614  ea^  =  -1° 

686.335  ea   =  -5° 


The  change  in  the  periods  is  clearly  stronger  for  a 
colatitudinal  initial  displacement  (c-,  ^   0,  see  Section  3.1) 
The  table  given  includes  nonlinear  effects,  but  for  small 
initial  displacement,  we  should  expect  that  most  of  the 
change  can  be  accounted  for  by  the  linear  theory. 
Figures  A. 71,  A. 72,  A. 73   show  some  examples  of  motions  for 
eBn  =  1°,   eB,  =  -2°,   ea,  =  -1°,  respectively. 
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(b)   8  =  49°,   c  =  .7.   The  configuration  is  0.6864° 

away  from  the  critical  value,  inside  the  linear  stable  zone. 

For  colatitudinal  initial  displacements.-  the  motions 
present  the  same  characteriztics  as  the  previous  configuration. 

When  the  initial  disturbance  is  in  the  longitudinal  direc- 
tion, the  situation  changes  radically.   The  amplitudes  of 
the  oscillations  of  the  colatudinal  components  are  about 
4  times  larger  than  the  amplitude   of  the  initial  displace- 
ment, while   D  (t) ,  k  =  1,...,4,  remains  bounded  by 
the  displacement. 

If  we  decrease  the  initial  displacement,  the  amplitudes 
decrease  in  approximately  the  same  amount,  and  hence  the 
motion  can  be  considered  as  a  stable  linear  motion.   In 
this  region,  nonlinear  effects  start  to  enter  into  the 
motion.   Some  runs  were  made  for   0  =  48°.   The  results  show 
that  the  nonlinear  effects  tend  to  disappear. 

In  Figure  A. 74  we  have  one  configuration  under  the 
longitudinal  perturbation  ea^  =  0.1°. 

The  change  in  the  angular  velocity,  presents  similar 
characteristics  as  in  the  previous  case  whenever  c,  7^  0 
or  c,  =  0  (see  3.1). 

Configuration  at  49°,   c  =  .7,   y  =  .15-10 

P 

Equilibrium  period   948.983  hours 

Nonlinear  Period  (Average)   Initial  Displacement. 

949.544   hours  ea^  =  0.1° 

948.918  ea^  =  -  0.5° 

1022.225  e^^   =        0.5° 

1051.566  e6j_  =  -  0.5° 
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(c)   i:  e  =50.41°,  c  =  .7,   The  configuration  is  ~  0.7° 
away  from  the  critical  value,  inside  the  exponentially 
unstable  zone. 

Figure  A. 75  shows  a  nonlinear  periodic  motion  with 
departure  away  from  the  equilibrium.   During  the  4000  hour 
integration,  the  motion  is  still  a  nonlinear  periodic 
motion  but  eventually  it  will  depart  away  from  the  neighbor- 
hood of  the  equilibrium  position. 

ii:  9  =  60°,  c  =  .5.   The  configuration  is  inside  the 
P 

exponentially  unstable  zone. 

Experiments  were  carried  out  only   with  initial 
displacements   eB,  ,  eot-,  between  -2.0°  and  2.0°. 

All  the  motions  obtained  were  nonlinear  periodic 
motions  with  departure  away  from  the  equilibrium  position. 
Figures  A, 76,  A. 77   show  the  resulting  motions,  under  a 
positive  eB-,  =  0.125°   and  negative  cB-^  =   -0.25°  colatitudi- 
nal  initial  displacement,  in  the  nonlinear  periodic  stage. 

From  the  overall  nonlinear  results  obtained,  we  can 
conclude  that  the  linear  analysis  gives  a  good  approximation 
to  the  motion.   Furthermore,  the  zones  of  exponential 
stability  and  instability  on  the  (9/C  )  plane   are  well 
defined  by  the  linear  theory  for  N  =  4,   in  the  sense  that 
whenever  the  linear  theory   establishes  an  exponentially 
stable  or  unstable  configuration,  it  is  confirmed  by  the 
nonlinear  calculations;  except  in  a  neighborhood  of  the 
boundary  of  the  linear  stable  zone.   This   neighborhood, 
for  N  =  4,  c   =  .7,   is  about  2°. 
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List  of  Figures 

Each  of  the  figures  (A. 71)  -  (A. 77)  consists  of  5 
subfigures  labeled   a,b,c,d,e.   The  subfigures  are  plots  of 
the  following  quantities : 


a:  6j_(t) 

b:  e2(t) 

c:  e3(t) 

d:  Q^U) 

e:  \(t)  ,  k 


=  1,2,3,4, 


For  example.  Figure  A. 72  c   represents  the  component  9^(t) 


Figure  Number 

Colatitude  of 
Circle  Vortices 

Coriolis  Para 

-meter,  c 
P 

Initial 
Displacement 

A. 71 

45° 

.7      1  e3j_  =   1.0° 

A. 72 

45° 

.7 

eB^  =  -2.0° 

A. 73 

45° 

.7 

ea^^  =  -1.0° 

A. 74 

49° 

.7 

ea^  =   0.1° 

A. 75 

50.41° 

.7 

e3^  =   0o5° 

A.  76 

60° 

.5 

e3^  =   0.125° 

A. 77 

60° 

.5 

£3-,^  =  -0.25° 
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(ii)  4  Vortices  with  Equal  Strength  u,    Equally  Spaced  on  a 
Circle  of  Colatitude  9,   Boundary  Condition   at  the 
Equator. 

Here  we  study  three  cases  for  the  coriolis  parameter 
c  =  .8.   According  to  the  linear  analysis  of  section  2, 
zero  is  an  eigenvalue  of  multiplicity  two.   Furthermore, 
for  6  <_  41.7070°,  the  configuration  is  exponentially  stable 
whereas  for  0  >_  41.7071°,  it  is  exponentially  unstable. 

Numerical  integrations  were  carried  out  for  values  of  9 
in  the  exponentially  stable  and  unstable  zone, 
(a)   9  =  45°,  c  =  .8.   The  configuration  is  exponentially 
stable. 

Numerical  integrations  were  conducted  for  initial 
displacements  eB-,  and  ea,   between  -20.0°  and  20.0°.   The 
results  show  that  all  the  motions  are  linear  stable  motions. 
For  all  the  experiments   D,  (t),  k  =  1,2,3,4,   is  bounded 
by  3/2  of  the  initial  amplitude  while  9,  (t) ,  k  =  1,2,3,4, 
by  the  initial  amplitude. 

When   the  initial  displacement  is  in  the  0  direction, 
the  circle  vortices  tend  to  oscillate  about  a  new  circle  of 
colatitude,  north  or  south  of  the  equilibrium  position, 
depending  on  the  sign  of  c6,. 

This  behavior,  as  in  the  case  with  no  boundary 
condition,  effects  the  angular  velocity  of  rotation. 

In  the  table  below  we  have  listed  the  nonlinear  periods 
of  rotation  for  different  initial  displacements. 
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Configuration  at  35°,   c  =  .8,   y  =  .15*10 
Equilibrium  period   327.573  hours 

(Average)  Initial  Displacement 


331.695  hours  eB,  =    .5 

456.822  eB^  =  15° 

323.508  eB,  =  -0.5° 

221.274  eB-[^  =  -15° 

327.525  ea^   =        .5° 

308.484  ea   =  15° 

327.480  ea^  =  -1.0° 

308.153  ea,  =  -15.° 


The  change  in  the  periods  of  rotation  is  clearly  stronger 
for  a  colatitudinal  displacement   (c,  ^   0,  Section3.1).  The 
table   given   includes  nonlinear  effects,  but  for  small  pertur- 
bation, most  of  the  change  can  be  accounted  for  by  the  linear 
theory. 

In  Figures  A. 78,  A. 79,  A. 80  we  have  some  results  on 
the  nonlinear  calculations  for  eB-i  =  0.5°,   eB-,  =  -15°, 
ea-,  =  -15°   respectively. 

(b)   e  =  41°,   c  =  .8.   The  configuration  is  0.7170°  away 
from  the  critical  value,  inside  the  exponentially  stable  zone. 

For  a  colatitudinal  displacement,  the  motions  present 
the  same  characteristics  as  the  previous  configuration. 
For  a  longitudinal  initial  displacement,  the  amplitudes 
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of  9j^(t)  ,  k  =  1,2,3,4,  are  about  four  times  the  initial  pertur- 
bation, but  as  ea,    decreases  in  magnitude,  the  amplitudes 
decrease  in  the  same  amount.    In  this  region,  nonlinear 
effects  enter  into  the  motion. 

o 

Some  experiments  were  conducted  for   9  =   40  .  They  show 
that  the  nonlinear  effects  that  produce  large  amplitudes  in 
the  colatitudinal  components  tend  to  disappear. 

(c)   0  =  50°,   c   =  .8.   The  configuration  is  in  the 

P 

exponentially  stable  range. 

The  experiments  were  done  using  perturbations  in  the 
positive  and  negative  colatitudinal  and  longitudinal  directions 
with  a  maximum  amplitude  of  0.5°. 

For  all  the  calculations,  the  motions  are  nonlinear 
periodic  motions  in  the  deparature  away  from  the  equilibrium 
position.   Figures  A. 81,  A. 82  and  A. 83   present  some  motions 
in  the  exponentially   unstable  range  for  e3-,  =  0.25°, 
ea,  =  0.25°,  eS-,  =  0.1°   respectively. 

From  the  results  of  the  numerical  integration  of  the 
nonlinear  equations  of  vortex  motion,  we  can  state  that  the 
linear  analysis  gives  a  good  approximation  to  the  motion. 
Furthermore,  the  zones  of  exponential  stability  and 
instability  on  the  (0,c  )  plane  are  well  defined  by  the 
linear  theory  for  N  =  4,  except  in  a  neighborhood  of  the 
boundary  of  the  exponentially   stable   range.   This  neighbor- 
hood is  about  1.5°  for  N  =  4  ,  c  ==  .8. 

P 
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List  of  Figures 

Each  of  the  figures  A. 78  -  A. 83   consists  of  5  sub- 
figures  labeled   a,b,c,d,e.   The  subfigures  are  plots 
of  the  following  quantities: 

a:  0^(t) 

b:  e2(t) 

c:  e3(t) 

d:  e4(t) 

e:  Dj^(t)  ,  k  =  1,2,3,4. 


For  example.  Figure  A. 72  c  represents  the  component 


e3(t) 


Figure 
Number 

Colatitude 
of  Vortices 

Coriolis  Parameter 

Initial 
Displacement 

A. 78 

35° 

.8 

e3-L  =    0.5° 

A. 79 

35° 

.8 

e&^  =   -15.0° 

A. 80 

35° 

.8 

ea^^  =  -15   ° 

A. 81 

42° 

.8 

e3j_  =     .25° 

A. 82 

42° 

.8 

ea^  =     .25° 

A. 83 

50° 

.8 

egj^  =    0.1  ° 
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